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Executive  Summary 

The  Final  Technical  Report  for  Grant  F49620-96- 1-0443  “Electron  Transfer  Chemistry  in  Optical 
Materials:  an  EPR  Investigation  of  Radiation-Induced  Defects  in  Chemically  Modified  Materials”  details 
the  construction,  and  performance-evaluation  of  a  94.9  GHz  EPR  (electron  paramagnetic  resonance) 
spectrometer,  operating  in  both  CW  (continuous  wave)  and  pulsed  modes.  The  first  application  of  this 
instrument  in  the  study  of  radiation  induced  defects  in  optical  materials  -  neutron-irradiated  a-Al203 
single  crystal  fiber,  and  y-irradiated  of  vitreous  silica  samples  -  are  also  described. 

A  major  finding  of  this  work  is  that,  owing  to  its  excellent  sensitivity  for  volume-limited  samples, 
high-frequency  EPR  instrumentation  makes  it  possible  to  characterize  paramagnetic  centers  in  optical 
fibers.  This  capability  is  illustrated  by  reported  spectroscopic,  line  shape,  and  analytical  studies  of  Cr3*, 
Fe3+,  and  color  centers  (electrons  trapped  at  anion  vacancies)  at  ppm  levels  in  an  a-Al203  fiber.  Special 
characteristics  of  high-magnetic  fields  for  improving  and  facilitating  EPR  spectroscopy  and  analysis 
(spin-counting)  are  detailed. 

An  emergent  feature  of  high-frequency  EPR  spectroscopy  is  the  apparent  ease  of  eliciting 
saturation  and  rapid  passage  effects.  Pioneering  efforts  in  the  systematic  characterization  and 
theoretical  simulation  of  these  effects  are  reported.  Preliminary  results  suggest  that  passage 
effects  can  be  harnessed  to  provide  enhanced  sensitivity  and  pronounced  relaxation-dependent 
selectivity  in  the  EPR  spectroscopy. 

The  personnel  that  participated  in  this  work  are: 

David  J.  Singel,  PI;  Vladimir  Krymov,  Visiting  Scientist  for  the  Physico-Technical  Institute  Donetsk; 
Dwight  A.  Schwartz,  postdoctoral  research  associate;  Eric  D.  Walter,  graduate  student  research  assistant; 
and  Sean  J.  Mcllwain,  undergraduate  research  aide. 

Publications  and  presentations  emanating  form  this  work  are  listed  below:. 

“High-Frequency  (94.9  GHz)  EPR  Spectroscopy  of  a  Spin  3/2  Center  in  a  Neutron-Irradiated  Sapphire 
Single  Crystal.”  D.A.  Schwartz,  Eric  D.  Walter,  Sean  J.  Mcllwain,  Vladimir  N.  Krymov,  and  David 
J.  Singel,  Appl.  Magn.  Res.  (1999)  16  223-236. 

“EPR  Studies  of  Optical  Materials:  Characterization  and  Control”  Invited  Lecture,  OpTeC  Conference  on 
Optical  Science  and  Laser  Technology,  Bozeman,  MT  September  1995. 

“Magnetic  Resonance  Spectroscopy  of  Optical  Materials”  Invited  Lecture,  OpTeC  Conference  on  Optical 
Science  and  Laser  Technology,  Bozeman,  MT.  August  1997. 

’’High  Field  EPR  Spectroscopy  of  Optical  Materials”  Invited  Lecture,  Optical  Science  and  Laser 
Technology  Conference,  Bozeman,  MT.  August  1998. 

“Experts  Workshop  on  High  Field  EPR  Spectroscopy”,  Workgroup  Leader,  Washington,  DC.  September 
1998. 

“High  Field  EPR  Spectroscopy  of  Optical  Materials  OpTeC  Conference  on  Optical  Science  and  Laser 
Technology,  Bozeman,  MT.  August  1999. 

“Passage  Effects  in  High-frequency  EPR  Spectroscopy  of  Optical  Materials”  OpTeC  Conference  on 
Optical  Science  and  Laser  Technology,  Bozeman,  MT.  August  2000. 

“Passage  Effects  in  High-frequency  EPR  Spectroscopy”  Invited  Lecture,  International  conference  on 
Magnetic  Resonance,  Florence,  Italy  2000. 
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1.  Introduction 


The  research  program  on  “Electron  Transfer  Chemistry  in  Optical  Materials  ”  at  Montana  State 
University  had  two  primary  objectives: 

•  Development  and  implementation  of  a  high-frequency  (94.9  GHz)  EPR  instrument; 

•  Application  of  high-frequency  EPR  spectroscopy  to  redox  active  centers  in  solid-state 
materials,  in  particular,  materials  of  interest  in  opto-electronics. 

This  report  provides  technical  details  on  the  instrument  that  has  been  constructed  at  Montana 
State  University,  and  results  obtained  in  the  initial  experimental  work  undertaken  with  the 
instrument. 

Since  its  earliest  days  of  application,  EPR  (electron  magnetic  resonance)  spectroscopy  has 
conventionally  been  performed  at  operating  frequencies  of  order  101  GHz.  Following  the 
pioneering  efforts  of  the  Lebedev  group  in  Moscow  [1],  there  is  currently  a  surge  of  interest  in 
the  development  and  application  of  EPR  spectroscopy  at  “high”  frequencies,  of  order  102  GHz. 
This  interest  is  more  developed  in  the  European  community  than  in  the  United  States.  There  is 
and  multinational  workgroup  on  high-frequency  EPR  in  Europe,  featuring  numerous  research 
groups  in  over  8  nations.  In  the  United  States,  by  contrast  progress  in  this  field  has  been  slow. 
The  system  implemented  at  Montana  State  University  is  only  high-frequency  instrument  west 
of  the  Mississippi  River,  and  one  of  only  a  few  instruments  in  the  country  capable  of  operation 
in  both  pulsed  and  CW  modes. 

The  leap  in  spectrometer  frequency  from  ~9.5  to  95  GHz,  with  its  corresponding  impact  on 
external  magnetic  field  strengths,  provides  a  powerful  means  to  alter  the  balance  between  the 
Zeeman  and  non-Zeeman  interactions  (electron-nuclear  hyperfine  and  electron-electron  fine 
structure  interactions)  of  a  spin  system,  and  thus,  to  alter  their  relative  importance  in  shaping  the 
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spectral  display.  An  important  goal  of  such  manipulations  is  to  forge  a  stronger  link  between  the 
features  of  a  spectrum  and  the  physical  quantities  to  be  measured.  For  example,  in  arguably  the 
most  fruitful  realm  of  application,  high  fields  are  used  to  magnify  the  spectral  effects  of  Zeeman 
anisotropy  in  free-radicals,  so  that  the  Zeeman  powder  pattern  is  prominently  featured  in  the 
spectral  display.  High-frequency  EPR  thus  promotes  direct,  precise  measurement  of  the 
principal  values  (g-values)  of  the  Zeeman  interaction,  and  thereby  facilitates  the  identification  of 
the  radical  and  the  characterization  of  its  environment  [1,2].  Conversely,  high  fields  have  also 
been  utilized  to  reduce  line-broadening  effects  in  the  central  transitions  of  Mn2+  ions,  thereby 
heightening  sensitivity  to  hyperfine  interactions  [3].  There  are  a  number  of  additional  important 
benefits  that  high-frequency  EPR  spectroscopy  can  afford.  Most  notably  is  the  promise  of  high 
sensitivity  in  volume  limited  samples.  This  feature  has  especially  important  ramifications  for  the 
study  optical  fiber  materials.  In  previous  EPR  studies  of  such  systems  investigators  studied 
“preforms”  -  the  bulk  material  from  which  fibers  are  drawn.  The  work  detailed  is  the  first 
demonstration  that  high-frequency  EPR  brings  high  sensitivity  to  the  study  of  paramagnetic 
centers  in  fibers,  including  both  metal  ion  inclusions  and  radiation  induced  “free-electron ” 
defects. 

An  final  feature  of  interest  is  the  utility  in  the  study  of  high-multiplicity  spin  systems  -  in 
particular,  systems  that  have  zero-field  splittings  which  are  comparable  to  the  microwave 
frequency,  v0,  of  conventional  spectrometers,  and  therefore  are  considerably  smaller  than  v0  in 

“high-frequency”  instruments.  The  work  presented  here  shows  that  high-frequency  EPR  leads 
to  a  simplification  of  the  otherwise  complicated  spectra.  Such  simplification  enables  improved 
quantification  both  of  the  number  of  paramagnetic  centers  and  the  parameters  that  determine  the 
spectrum. 
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2.  Instrumentation 


The  high-frequency  EPR  instrument  that  was  constructed  at  Montana  State  University  entails: 

•  an  Oxford  Instruments  Teslatron  split-pair  superconducting  magnet  system  with 
“intelligent”  power  supply  (IPS-120-10)  for  operation  from  0-3.7  T; 

•  a  94.9  GHz  microwave  bridge  obtained  from  collaborators  at  the  Physico-Technical 
Institute  of  Donetsk,  Ukraine; 

•  a  TEoi  1  cavity  assembly  made  by  collaborators  at  the  Physico-Technical  Institute  of 
Donetsk,  Ukraine; 

•  a  probe  -  constructed  locally  -  that  links  the  bridge  and  cavity,  and  provides  for  cavity¬ 
tuning,  microwave  coupling,  sample  rotation,  and  field  modulation.  The  installation  of  the 
bridge,  and  the  design  and  implementation  of  the  probe  was  done  by  Vladimir  Krymov  of 
Donetsk,  while  at  Montana  State  University. 

The  magnet  system  allows  for  computer-controlled  sweeping  of  the  field  exerted  by  the  split- 
pair,  over  an  operating  range  of  0  -  3.7  T,  with  a  maximum  rate  of  0.33  T/min  and  a  resolution  of 
3.3  pT.  Further  details  of  the  system  are  available  from  Oxford  Instruments.  The  fixed  current 
leads  of  the  magnet  are  optimized  for  low  helium  consumption  while  sweeping  the  field  through 
its  full  range.  For  many  spin  systems,  in  particular  those  with  appreciable  fine- structure 
interactions,  high-frequency  EPR  spectroscopy  requires  wide-ranging  field  sweeps.  The 
requirement  argues  against  magnet  systems  in  which  current  is  swept  in  an  auxiliary,  rather  than 
the  main,  magnet  coil.  While  such  dual  coil  systems  are  thought  to  enjoy  a  compelling 
advantage  in  substantially  reduced  helium  consumption,  the  magnet  system  implemented  in  this 
work,  with  its  optimized  current  lead  configuration,  gives  satisfactory  performance:  less  than  100 
liters  were  required  during  the  first  month  of  continuous  operation.  The  magnetic  field  was 
calibrated,  for  these  measurements,  by  the  method  of  Burghaus  and  Mobius  [13].  The  calibrated 


3 


value  of  the  field  consistently  lay  within  9  G  of  the  nominal  value  on  the  Oxford  Teslatron 
display. 

The  bridge  features  a  highly  stabilized  oscillator  (STAMO)  with  a  fixed  frequency  of  94.9  GHz; 
operation  entails  both  pulsed  and  CW  (continuous-wave)  modes.  The  maximum  microwave 
field  in  the  cavity  is  ~3.6  G  with  incident  power  of  250  mW.  When  operating  in  CW  mode,  the 
maximum  power  is  limited  by  fixed  attenuation  to  5.5  mW.  A  heterodyne  detection  scheme  is 
employed  in  the  bridge,  with  dual  outputs  in  phase-quadrature;  phase  settings  are  continuously 
variable  by  manual  adjustment.  Further  details  of  the  bridge  design,  including:  principles  of 
operation,  block  diagram  of  the  microwave  assembly,  lists  of  the  electronic  components,  and 
performance  specifications  are  included  in  Appendix  I.  Field  modulation  at  frequencies  up  to 
150  kHz,  and  maximum  amplitude  of  15  G  can  be  applied;  two  computer-interfaced  lock-in 
amplifiers  (EG&G  PAR  models  5207  and  5208)  are  used  to  detect  the  modulated  signal. 

The  cavity  assembly  is  illustrated  on  the  following  insert  pages.  The  first  shows  a  schematic  of 
the  assembly,  which  features  a  (slotted)  TE0n  right  circular  cylindrical  side-wall  coupled  to  the 
waveguide  transmission  line  through  an  iris.  As  indicated  in  the  illustrations,  the  waveguide  is 
part  of  an  assembly  that  remains  “fixed”  in  the  magnet,  while  the  cavity  is  part  of  an  assembly 
that  is  removable,  thus  facilitating  sample  changing.  Both  parts  are  designed  to  be 
accommodated  within  an  Oxford  Instruments  CF935G  dynamic  continuous  flow  cryostat 
(operating  at  temperatures  as  low  as  3.8  K).  The  cavity  is  spring  mounted  on  the  probe  to  insure 
an  appropriate  fit  of  the  coupling  components.  Coupling  is  adjusted  varying  the  by  orientation  of 
the  cylinder  axis  of  the  microwave  cavity.  This  movement  is  accomplished  by  rotating  the  entire 
removable  assembly  by  means  of  the  arm  at  the  top  of  the  probe;  the  alignment  of  this  axis  can 
be  varied  through  ±7t/8  from  its  nominal  orientation,  perpendicular  to  the  external  field.  The 
resonant  frequency  for  the  TEoi  i  model  of  the  cavity  is  made  to  match  the  94.9  GHz  oscillator 
frequency  by  remote  movement  (knurled  knob  control)  of  one  of  the  cavity  endwalls.  The 


4 


Bo 


Schematic  of  the  microwave  cavity  assembly 

The  microwaves  are  coupled  from  the  W-band  transmission  line  (C)  through  the  iris  coupling  hole  into  the 
right  circular  cylindrical  resonator,  which  is  slotted  to  enable  penetration  of  the  of  the  oscillating  rf  fields 
(depicted,  A)  or  audio  fields  for  field  modulation  (not  depicted).  The  sample  is  introduced  through  holes 
centered  in  the  plungers  (B),  whose  lateral  position  is  adjust  to  match  the  cavity  resonance  and  fixed 
microwave  frequency. 


Views  of  the  exterior,  fixed  assembly  portion  of  the  EPR  probe 

(A)  Bottom  end  showing  the  waveguide  transition  (f)  and  bend  to  reach  the  microwave  cavity  (b).  A 
temperture  sensing  element  is  mounted  at  (d). 

(B)  Top  end  showing  connection  of  waveguide  to  the  microwave  bridge  (a),  clearance  for  electrical 
connections  (c)  and  arm  (e)  for  adjusting  the  orientation  of  the  cavity  and  the  microwave  coupling. 
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Views  of  interior,  removable  portion  of  the  EPR  probe 

(A,  B)  Views  showing  the  cavity  assembly  (B)  and  its  extension  (A)  toward  the  top  of  the  probe.  The 
internal  rod  (b),  with  a  cone  at  its  bottom,  pushes  against  lever  (d)  thus  adjusting  the  position  of  the  cavity 
plunger.  The  clamped  (h)  sample  holder  (g)  is  rotated  by  a  shaft  (f)  and  gear  assembly. 

(C)  Views  of  the  top  of  the  probe  showing  adjustment  knobs  for  control  of  cavity  tining  (a)  and  sample 
rotation  (e)  with  respective  control  rods  (b  and  f).  Inputs  and  cables  (i)  for  rf  and  audio  fields  are  shown. 


position  of  the  other  wall  is  varied  by  the  introduction  of  spacers  to  keep  the  coupling  hole 
approximately  in  the  center  of  the  cavity.  This  adaptation  allows  for  the  introduction  of  a 
rotating-shaft  and  gear  assembly  for  sample  rotation  about  the  cylinder  axis  of  the  cavity,  as 
indicated  in  the  figures.  Cables  and  wiring  for  magnetic  field  modulation  (and  ultimately  for 
ENDOR  (electron  nuclear  double  resonance),  as  depicted)  and  for  temperature  sensing  are  also 
illustrated.  The  basic  features  this  system,  in  particular  the  remote  adjustable  tuning  and 
coupling  schemes  were  developed  over  years  of  collaboration  between  Lebedev  and  KrymovJ 

The  spectrometer  is  operated  with  software  that  was  written  under  the  Lab  View  programming 
package  at  Montana  State  University.  This  software  systems  provides  for  control:  of  the 
magnetic  field  stationary  setting,  or  field  sweep  between  specified  values  at  a  specified  rate;  the 
magnetic  field  modulation  frequency  and  amplitude,  and  the  settings  (phase,  gain,  harmonic, 
offset,  time-constant)  of  lock-in  amplifiers,  when  the  system  is  operating  in  CW  mode.  In  pulsed 
mode  the  master  clock,  pulse  forming  units  (Stanford  Research  Systems,  DG535),  boxcar 
integrator  (SRS  250/280),  and  D/A  conversion  system  (SRS  245)  are  software  controlled.  The 
software  system  provides  a  convenient  user  interface  for  adjustment  of  the  various  control  and 
acquisition  settings,  and  for  display  and  manipulation  of  the  output  EPR  spectra.  This  computer 
display  is  illustrated  on  the  accompanying  insert  pages.  The  LabView  “code”  is  reproduced  in 
Appendix  II. 

3.  Experimental  and  Results 

3.1.  Neutron-irradiated  sapphire  single-crystal  fiber 

Interest  in  this  system  -  as  part  of  a  broader  study  of  electron-transfer  chemistry  in  optical 
materials  -  is  sparked  by  the  photobleaching  and  thermoresistant  spectral  hole-burning 
phenomena  that  it  exhibits  [4,5].  These  phenomena  have  attracted  interest  to  the  material  as  a 
candidate  for  application  in  optical  data-storage  technologies  based  on  spectral  hole-burning  [6]. 
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e/Time 


Labview  Routine  for  EPR  Control  and  Data  Acquisition:  Experiment  Set-Up  Screen 


Labview  Routine  for  EPR  Control  and  Data  Acquisition:  Experiment  In-Progress  Screen 


Neutron-irradiation  of  ostensibly  white  sapphire,  followed  by  thermal  annealing,  generates  a 
series  of  new  optical  absorption  bands  [7].  Laser  irradiation  in  these  bands  can  produce  a 
frequency-selective  bleaching  of  the  irradiated  band,  accompanied  by  an  intensity  increase  in 
other  bands  [8],  This  process  has  been  modeled  as  photo-induced  electron  transfer  among  mono- 
and  di-anion  vacancies  with  one  or  two  electrons  trapped  at  each  vacancy  -  i.e.,  among  F  and  F2 
centers  having  various  “oxidation-states”  (0,  +1,  +2)  and  spin  multiplicities  (S=0,  1/2,  1,  3/2,  2), 
as  tabulated  below  [7,8],  Compelling  structural  evidence  for  this  model,  however,  is  lacking. 
Since  EPR  spectroscopy  is  well-established  as  a  choice  method  for  identifying,  characterizing, 
and  quantifying  paramagnetic  color  centers  -  and  thus  for  providing  crucial  information  required 
for  understanding  the  performance,  and  for  optimizing  the  preparation  of  such  materials. 
Motivated  both  by  the  photo-chemical  properties  of  neutron-irradiated,  thermally  annealed 
sapphire,  and  its  potential  applications,  we  launched  an  EPR  investigation  aimed  at  determining 
the  nature  of  the  paramagnetic  centers  in  this  material. 


Model  for  Induced  Centers: 


Oxygen  Vacancies 

Trapped  Electrons 

S 

^ma\  (nitl) 

F+ 

1 

1 

1/2 

205 

F 

1 

2 

0 

258 

f2++ 

2 

2 

1 

450 

f2+ 

3 

3/2 

355 

f2 

2 

4 

0 

300 

Samples  were  obtained  from  Professor  Aleksander  Rebane  at  Montana  State  University. 

Crystals  were  prepared  by  Professor  William  Yen  (University  of  Georgia)  by  a  laser-heated 
pedestal  growth  method,  with  starting  materials  of  99.999%  purity  [12].  This  method  enjoys  a 
smallness-of-scale  that  facilitates  rapid  preparation  of  a  wide  variety  of  crystal  samples.  Crystals 
fabricated  by  this  method,  however,  are  generally  small:  our  sapphire  samples  have  the  shape  of 
a  right  circular  cylinder,  0.5  mm  in  diameter  and  10  mm  in  length.  Note  that  these  fiber-like 
samples  have  a  volume  that  is  a  factor  of  10r  smaller  than  sample  volumes  typically  enjoyed  in 
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X-band  experiments.  Crystals  were  neutron-irradiated  with  a  dose  of  ~1019/cm2  and  annealed  at 
-600  K  by  Professor  Ilmo  Sildos  (University  of  Tartu,  Estonia)  [4,5]. 

3.1.1.  Spectroscopy 

Initial  results  of  the  high-frequency  EPR  investigation  on  this  material  have  been  published  in 
Ref.  27.  Three  predominant  sets  of  EPR  signals  (evident  in  Figure  1),  were  readily  identified  as 
Cr3+  (S=3/2)  and  a  pair  of  crystallographically  equivalent  Fe3+  (S=5/2)  impurity  ions.  High- 
frequency  EPR  instrumentation  played  a  key  role  in  facilitating  the  identification  of  the  high  spin 
multiplicity  Cr3+and  Fe3+  centers.  Numerous  weaker  EPR  signals  are  also  apparent  in  spectra 
displayed  with  expanded  vertical  scales  (see  Figure  2).  These  spectra  presumably  belong  to  the 
neutron  induced  color  centers,  and  have  subsequently  been  more  carefully  investigated  (as 
exemplified  in  Figure  3);  the  spectra,  as  detailed  below,  are  complex;  a  full  analysis  is  underway 
now  that  the  major  signals  have  been  identified. 

In  Ref.  27,  the  following  EPR  parameter  values  were  published: 

Cr3+  g,,  =  1.9812±0.0001,  |D|  =  5.738±0.003  GHz,  g±  =  1.9814±0.0001; 

Fe3+  g  =  2.0034±0.0001  and  |D|  =  5.033±0.002  GHz,  and  |a  -  F|  =  0.986±0.006  GHz. 

The  estimated  errors  are  evaluated  by  propagation  of  an  estimated  error  in  Bres  of  1 .7  G,  which 

derives  from  errors  in  observation  and  calibration  [13]. 

Comparison  of  the  amplitudes  of  the  Fe3+  or  Cr3+  EPR  signals  (adjusted  for  their  different 
values  of  S  and  S±)  with  those  measured  under  equivalent  conditions  for  an  EPR  dosimetry 

standard  (E’  centers  in  “wet”  suprasil  silica  exposed  to  a  dose  of  3  Mrad  60Co  y-irradiation)^, 
gives  an  experimental  estimate  of  1016  ions/cm3  for  both  Fe3+  and  Cr3+.  This  concentration 
corresponds  to  a  substation  ratio  of  10'6,  a  value  similar  to  the  purity  of  starting  material. 
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Figure  1.  94.9  GHz  EPR  spectrum  of  a  neutron-irradiated  sapphire  single  crystal  recorded  at  room  temperature  and 
with  the  external  field  making  an  angle,  ©=0.18rt ,  with  respect  to  the  crystal  c  axis.  The  microwave  power  was  3 
pW.  The  field  was  swept  at  0.33  T/min  and  modulated  at  75  kHz  with  an  amplitude  of  10  G.  Lock-in  detection 
was  employed  with  a  time-constant  of  100  msec. 


Figure  2.  94.9  GHz  EPR  spectrum  of  a  neutron-irradiated  sapphire  single  crystal  recorded  at  room  temperature  and 
with  the  external  field  at  ©=0. 1  8ti.  The  microwave  power  was  5  nW.  The  field  was  swept  at  100  G/min  and 
modulated  at  75  kHz  with  an  amplitude  of  1.5  G.  Lock-in  detection  was  employed  with  a  time-constant  of  100 
msec. 
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Figure  3.  Detail  of  94.9  GHz  EPR  spectra  of  a  neutron-irradiated  sapphire  single  crystal  recorded  at  room 
temperature  and  with  the  external  field  making  the  designated  angles  (degrees)  with  respect  tot  he  crystal  c  axis. 

The  microwave  power  was  3  nW.  The  field  was  swept  at  lOOG/min  and  modulated  at  100  kHz  with  an  amplitude  of 
1  G.  Lock-in  detection  was  employed  with  a  time-constant  of  30  msec. 

For  Cr3+,  Zverev  and  Prokhorov  previously  reported  g„  =  1.982±0.002,  g±  =  1.979±0.009,  ana 
D  =  -  5.73±0.03  GHz  [14];  Kornienko  and  Prokhorov  [15]  reported  for  Fe3+  the  parameter 
values  - 

g  (isotropic)  =  2.0030±0.0006,  D  =  +5.034±0.001  GHz,  and  |a  -  F|  =  0.990±0.001  GHz.  A 
remark  is  in  order  concerning  the  Cr3+  parameters  values.  Manenkov  and  Prokhorov  [22] 
originally  reported  the  g-values,  g„  =  1.984,  g±  =  1.9867  based  on  work  at  8.96  and  1 1.97  GHz. 

The  values  cited  above  were  obtained  later  in  “more  complete”  work  at  37.86  GHz  [14].  In  later 
work  by  several  other  groups  [23-25],  values  closer  to  the  original  ones  [22]  appear.  The 
suggestion  was  made  in  Ref.  27,  that  the  conventional  spin-Hamiltonian  used  to  describe  the 
Cr3+  might  not  be  entirely  adequate.  In  particular,  “high-spin  Zeeman”  terms  linear  in  field  and 
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cubic  in  spin  angular  momentum  -  which  are  admissible  in  the  most  general  spin-Hamiltonian 
for  an  S=3/2  system,  but  are  conventionally  ignored  -  had  been  suspected  as  the  origin  of  this 
discrepancy.  A  detailed  analysis  of  the  Cr^+  spectra  was  subsequently  undertaken  by  us  to 
address  this  issue;  elements  of  this  work  are  given  in  Appendix  III.  The  major  conclusion  is  that 
the  high-field  Zeeman  terms  are  negligible,  and  the  discrepancy  simply  emanates  from  some 
author’s  repeated  quotation  of  Prokhorov’s  original,  less  accurate  results. 

The  presence  of  transition  metal  ions  in  this  material,  at  concentrations  that  may  be  simlar  to  that 
of  the  color  centers,  raises  additional,  fundamental  issues.  In  their  pioneering  study  of  the  optical 
absorption  spectra  of  neutron  damaged  sapphire  [7],  Evans  and  Stapelbroek  found  transition 
metal  ion  concentrations  to  be  less  than  ~5  ppm.  They  provided  detailed  arguments  that  the 
absorption  spectra  were  not  assignable  to  substitutional  metal  ion  impurities,  but  neither  they  nor 
subsequent  investigators  fully  considered  the  effect  that  their  presence  might  have  on  the 
production  or  aggregation  of  lattice  defect  centers.  Redox  active  ions  should  be  expected  to 
impact  the  solid-state,  electron-transfer  chemistry  that  occurs  in  the  preparation  (irradiation  and 
thermal  annealing)  and  the  utilization  (optical  photo-bleaching,  hole-burning)  of  this  material. 

For  example,  a  recent  study  of  the  effects  of  neutron-radiation  on  metal  ion  doped  sapphires 
clearly  shows  that  the  presence  of  the  ions  alters  the  distribution  among  radiation-induced 
products;  in  particular,  new  products  involving  F  centers  associated  with  metal  ions  are  reported 
[19].  Dose-response  effects  are  modulated  by  the  radiation-induced  valence  changes  in  Cr  and 
Fe  centers  -  an  effect  also  reported  in  y-irradiation  studies  [20],  Such  effects,  while  not  well- 
understood  systematically,  have  nevertheless  been  exploited  empirically  to  impart  radiation 
hardness  to  optical  materials  [21].  For  photo-ionization  and  electron-transfer  processes,  the 
potential  effects  of  redox  active  are  even  more  obvious.  Studies  of  the  optical  properties  of 
neutron-irradiated  sapphire  must  include  an  assessment  of  transition  metal  content,  and  must  be 
designed  to  control  for,  and  to  exploit,  their  effects. 
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EPR  Signal 


3.1.2.  Line  shapes  of  the  Chromium(III)  centers. 

An  analysis  of  the  line  shapes  of  the  Cr3+  EPR  signals  discussed  above,  demonstrated  (see 
Figure  4)  that  the  line  shapes  adhere  very  closely  to  (derivative)  Lorentzian  form  (width  of  ~30 
MHz)  and  deviate  noticeably  from  Gaussian  shapes.  This  result  seemed  surprising.  Inasmuch  as 

the  measured  Cr3+  spin  concentration  in  the  sample  is  1016/cm3,  the  interaction  with 
neighboring  aluminum  nuclear  moments  would  be  expected  to  be  the  dominant  line  broadening 
effect.  This  interaction,  however,  would  certainly  lead  to  an  inhomogeneous  Gaussian 

broadening  of  the  lines.  Cr3+'Cr3+  dipolar  interactions  could  lead  to  spin-exchange  line 
broadening  effects,  but  concentrations  of  102%m3,  would  be  required  to  produce  the  observed 
line  widths. 


Magnetic  Field  (Tesla) 


Magnetic  Field  (Tesla) 


Figure  4.  Line  shape  analysis  of  Cr3+  EPR  spectra.  The  same  experimental  spectrum  (Ms=±l/2  transition  with  the 
field  aligned  along  the  crystal  c  axis.)  is  depicted  in  both  the  right  and  left  panels.  The  least  squares  best  fit 
derivative  Gaussian  (left  panel)  and  Lorentzian  (right  panel)  are  superimposed  (dashed  lines)  on  the  spectra.  The 
Lorentzian  shape  gives  an  excellent  representation  of  the  data,  with  a  line  width  of  19  G). 


These  results  can  be  explained  by  assuming  that  the  distribution  of  the  impurity  centers  is  not 
uniform  within  the  lattice;  the  centers  are  segregated  from  the  bulk  material.  The  global 
concentration  (measured  by  spin  counting)  could  thus  remain  small,  while  the  local 
concentration  (reflected  in  the  line  shape)  could  be  very  large,  as  dictated  by  the  degree  of 
segregation.  To  test  this  idea,  surface  specific  elemental  analyses  -  SIMS  (secondary  ion  mass 
spectrometry)  and  XPS  (x-ray  photoelectron  spectroscopy  -  were  carried  out  at  the  ICAL  facility 
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of  Montana  State  University.  The  results  of  these  and  the  EPR  measurements  are  given  in  the 
table  below. 


Analysis  Method 

EPR  intensity 

EPR  line  width 

SIMS 

XPS 

Range  of  Sample 
Probed 

Global 

Near-neighbor 

Spins 

Depth  of -10  A 

Depth  of  -500  A 

Cr  Concentration 

1016/cm3 

6  x  1 020/cm3 

4  x  1021 

1  x  1020 

The  SIMS  and  XPS  support  the  idea  of  that  the  chromium  ions  are  substantially  segregated  and 
further  reveal  that  they  are  congregated  near  the  very  surface  of  the  fiber.  It  remains  to  be 
determined  whether  this  clustering  occurs  during  the  laser  pedestal  crystallization  process,  or 
during  the  subsequent  irradiation  and  annealing  of  the  particular  sample  studied.  It  has 
previously  been  noted,  however,  that  the  optical  T2’s  of  materials  grown  by  laser  pedestal 
methods  tend  to  be  markedly  smaller  than  those  of  materials  made  by  slower  methods  of  growth 
methods.  It  is  therefore  tempting  to  suggest  that  the  segregation  occurs  during  the 
crystallization  process:  this  partitioning  would  explain  both  the  anomalous  EPR  and  optical 
homogeneous  line  widths. 


3.2.  y-Irradiated  Vitreous  Fused  Silica 

The  paramagnetic  defects,  generically  known  ac  E’ -centers,  induced  by  y-irradiation  of  fused 
silica  glassy  samples  have  been  studied  in  detail  for  many  years  [29].  A  study  of  a  sample  of 
such  material  was  undertaken  with  one  major  goal  in  mind:  to  establish  the  performance  of  the 
instrument  built  in  the  time-domain ,  spin-echo  mode  of  operation.  Further  study  of  this  material 
was  warranted  by  the  opportunity  to  utilize  it  in  a  program  of  high-frequency  EPR  dosimetry, 
with  the  goal  of  assaying  novel,  biological  radiation  shielding  material,  as  reported  in  the  NASA 
innovation  magazine  [30].  The  CW  EPR  radiation  dosimetry  studies  led  to  an  appreciation  of 
prevalence  of  “rapid  passage”  effects  in  high-frequency  EPR  spectroscopy,  and  the  emergence  of 
a  new  line  of  investigation  aimed  at  elucidating  and  harnessing  these  effects. 
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Magnetic  Field  (T) 


Figure  5.  Echo-detected  EPR  spectrum  of  60Co  y-irradiated  suprasil  fused  silica  glass.  Experimental  details 
are  given  in  the  text.  Absent  field-modulation  and  lock-in  detection  the  integral  form  of  the  spectrum  is 
registered.  «• 


3.2.1.  Benchmarking  of  spin  echo  capabilities. 

E’-centers  in  y-irradiated  fused  silica  glass  are  a  standard  test  sample  for  evaluating  pulsed  EPR 
instrumentation.  An  echo-detected  EPR  spectrum  of  such  a  sample  is  shown  in  Figure  5.  The 
sample  is  a  0.35  mm  diameter  rod  of  suprasil  quartz  subjected  to  100  Mrad  of  60Co  y-irradiation. 
The  echo-detected  method  employed  involves  application  of  two,  94.9  GHz  microwave  pulses 
with  peak  power  -230  mW,  duration  36  nsec  separated  by  a  delay  of  200  nsec  and  detection  on 
the  resulting  echo  at  a  subsequent  delay  of  -  200  nsec  with  a  boxcar  integrator  whose  gate  width 
was  -50  nsec.  The  echo  cycle  was  repeated  at  a  rate  of  1  kHz.  The  sample  was  at  room 
temperature  in  these  experiments.  For  echo-detected  EPR  spectroscopy,  the  sequence  is 
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continually  repeated  as  the  magnetic  field  is  scanned.  An  exponential  average  of  the  prior  300 
echoes  was  acquired  at  each  field  step  to  compose  the  spectrum.  Single-shot  signal-top-noise 
ratios  were  102  under  these  conditions.  This  result,  is  similar  to  performance  standards  at  lower 
microwave  frequencies  (e.g.,  X-band  (-10  GHz) ,  provided  the  sample  is  increased  to  its 
optimum  size  at  the  lower  frequency  (101  greater  diameter  and  101  greater  length).  This  result 
thus  indicates  a  improvement  in  sensitivity  of  as  much  as  1 03  at  95  GHz,  for  volume  limited 
samples. 


3. 2. 2  Saturation  and  Passage  effects  in  high-frequency  EPR 

CW  EPR  spectra  of  the  same  fused  silica  sample  are  shown  in  Figure  6.  The  two  phase- 
quadrature  microwave  signal  channels  are  displayed  (absorption  signal,  solid  line;  negative 
dispersion  signal  dashed  line)  together  with  the  vector  magnitude  response  (dotted  line).  Under 


Figure  6.  Absorption,  Dispersion,  and  Magnitude  EPR  spectra  of  E’  center  in  irradiated  quartz. 
Microwave  frequency  94.9  GHz,  power  55  nW  (-50.6  dB  attenuator  setting),  scan-rate  200  G/min, 
time  constant  0.1  sec.  modulation  frequency  75  kHz,  modulation  amplitude  2.5  G. 
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EPR  Signals 


non-saturating  conditions,  the  absorption  and  dispersion  signals  comprise  a  simple  Hilbert 
transform  pair  -  as  is  neatly  illustrated  in  Figure  6.  If  the  microwave  power  is  increased 
sufficiently,  this  simple  relationship  is  broken.  The  absorption  signal  amplitude  reaches  a 
stationary  level  (saturation);  moreover  an  apparent  partial  integration  of  the  derivative  (field- 
modulated  and  lock-in  detected)  CW  EPR  signals  can  occur,  as  vividly  illustrated  in  Figure  7 


Magnetic  Field  (T) 

Figure  7.  Power  dependence  of  phase  quadrature  EPR  signals  of  the  irradiated  quartz  sample  (Figs.  5 
and  6).  The  detector  phase  is  adjusted  such  that  at  non-saturating  powers  (bottom  spectrum)  the  solid  line 
is  the  absorption  and  the  dashed  line  is  the  dispersion  spectrum  (Fig.  6).  Acquisition  parameters  are  the 
same  as  in  Fig.  6,  apart  from  microwave  power,  which  varies  top-to-bottom  )  as:  -1 1.7,  -28.9,  -39.0,  - 
42.8,  -45.9,  -48.5,  and  -50.6  dB  from  full  power  (5  mW). 
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[31].  (Only  the  passage  effects  are  illustrated;  the  signals  are  “peak  normalized”  thus  relative 
amplitude  information  is  not  exhibited).  An  intriguing  aspect  of  high-frequency  EPR 
spectroscopy  is  the  apparent  ease  of  reaching  saturation  in  the  study  of  dilute  paramagnetic 
centers  in  solids.  To  operate  in  the  linear,  non-saturating  it  has  been  necessary  to  operate  the 
instrument  built  MSU  at  incident  microwave  powers  of  ~nW,  far  from  its  5  mW  peak  CW 
power.  Farrar  et  al.  [32]  in  a  140  GHz  study  of  a  tyrosyl  radical  in  a  photosynthetic  reaction 
centers,  found  that  the  dispersion  signal,  in-phase  with  the  field  modulation,  obtained  under  rapid 
passage  conditions,  enabled  an  crucial  enhancement  of  signal-to-noise  ratio  over  the 
conventional  absorption  signal,  and  also  provided  better  resolution  (less  power-broadening)  of 
hyperfine  structure.  Lebedev  and  co-workers  [33]  have  also  published  140  GHz  EPR  spectra 
obtained  as  conventional  absorption  signals,  as  well  as  both  dispersion  signals  in  phase  and  in 
phase-quadrature  with  the  field  modulation,  and  by  analysis  of  these  signals  have  obtained 
estimates  of  characteristic  relaxation  times. 

Studies  of  this  behavior  have  a  number  of  antecedents  in  solid-state  EPR  at  conventional 
frequencies  (not  to  mention  the  entire  field  of  saturation-transfer  EPR  in  fluids  [34]).  In 
pioneering  work,  H>  ie  used  rapid  passage  methods  discern  the  presence  distinct  color  centers 
with  overlapping  spectra  in  irradiated  LiF  [35].  Very  recently  Godbole  and  Sastry  [36]  revisited 
this  idea,  and  demonstrated  that  certain  spectral  displays  can  be  used  to  “edit”  complex  EPR 
spectra  involving  multiple  paramagnetic  centers  with  distinct  relaxation  properties:  specifically, 
in  a  crystal  of  CaF2  co-doped  with  U3+  and  Gd3+  they  established  conditions  under  which  an 
overlapping  U3+  spectrum  was  edited  from  the  spectral  display,  leaving  a  “clear  Gd3+  spectrum”. 
Mailer  and  Taylor  [37]  found  a  signal  enhancement  of  at  least  500  in  ferricytochrome  c  by 
detecting  rapid  passage  EPR  signals  in  the  dispersion  mode.  Ammerlaan  and  van  der  Wiel[38] 
measured  longitudinal  relaxation  times  of  vacancies  in  silicon  by  analysis  of  the  temperature 
dependence  of  dispersion  EPR  lineshapes  obtained  under  passage  conditions.  Such  examples 
make  it  clear  that  there  are  significant  advantages  to  be  gained  in  signal  strength,  spectral 
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selectivity,  sensitization  to  dynamics  that  can  be  accessed  by  departing  from  the  conventional 
practice.  At  the  same  time  it  must  be  admitted  that,  in  contrast  to  the  study  of  dynamics  in  fluids 
by  saturation-transfer  EPR,  the  application  of  these  potential  tools  in  the  EPR  of  solids  has  been 
limited. 

In  part,  this  limitation  is  arguably  related  io  instrumentation.  Good  signal-to-noise  ratios  in 
dispersion  EPR  are  hampered  by  demodulation  of  source  noise.  This  problem  has  been 
mitigated  by  the  use  of  a  bimodal  resonator[39]  and  high-bandwidth,  loop-gap  resonators[40]. 
High-frequency  resonators  with  quality  factors  of -5000  enjoy  a  similar  bandwidth  to  low  Q 
loop-gap  resonators,  and  are  thus  likewise  well-suited  to  dispersion  detection.  A  second  aspect 
can  be  linked  to  the  complexity  of  the  underlying  theory.  The  Bloch  equations  augmented  with 
sinusoidal  modulation  of  the  magnetic  field  are  not,  in  general,  analytically  solvable.  Mortis  [31] 
and  Weger  [32]  derived  certain  expressions  that  are  applicable  in  certain  limiting  cases 
(although  not  always  experimentally  useful  ones).  Their  work  (reprised  also  by  Ammerlaan  and 
van  der  Weil  [38]  elucidated  much  of  the  basic  phenomena,  in  particular  the  nature  of  the 
absorption  and  dispersion  signals  obtained  in  phase  and  in  phase  quadrature  with  the  field 
modulation.  Cullis  used  a  spin  temperature  approach  in  treating  this  problem  [43].  He 
specifically  considered  issues  of  spectral  line  width,  and  power-broadening;  moreover,  on  the 
basis  of  his  analysis,  he  suggested  new  techniques  for  measuring  relaxation  properties.  Solutions 
to  this  problem  have  been  advanced  —  especially  in  the  realms  of  coherent  optics  and  saturation- 
transfer  EPR  -  involving  Fourier  or  Floquet  [44]  expansions,  followed  by  evaluation  of  the 
expansion  coefficients  by  continued  fraction,  recursive,  and  other  mathematical  approaches 
[34?45_46].  Grivet  [47]  carried  out  numerical  solutions  of  these  augmented  Bloch  equations  by 
direct  numerical  integration  of  the  equations  using  the  fourth-order  Runge-Kutta  algorithm  [48]. 
While  computational  solutions  such  as  these  do  not  directly  advance  our  understanding,  they  do 
provide  a  valuable  tool  to  simulate  the  behavior,  and  thereby  providing  conceptual  support  to  the 
experimental  work,  and  enabling  its  extension  to  situations  not  yet  encountered  experimentally. 
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Rapid  passage  techniques  in  solid-state  EPR  spectroscopy  thus  appears  to  have  the  potential  to 
yield  large  enhancements  of  signal  amplitudes,  and  to  provide  selectivity  that  can  be  used  for 
spectral-editing.  Moreover,  passage  phenomena  are  particularly  prevalent  in  high-frequency 
EPR.  Accordingly,  a  pilot  program  combining  high-frequency  EPR  spectroscopy  of  test-bed 
materials  and  numerical  integration  of  the  underlying  equations  of  motion  was  initiated  to 
explore  the  feasibility  of  realizing  this  potential.  The  test-bed  centers  studied  the  E’  center 
formed  in  irradiated  vitreous  silica,  and  the  free  radical  BDPA  (a,y-bisdiphenylene*P- 
phenylallyl)  incorporated  into  a  matrix  of  atactic  polystyrene  (PS)  -  the  same  sample  employed 
by  Farrar  et  al.[32] 


In  the  numerical  work,  the  simple  approach  of  Grivet  [47]  was  followed,  namely  numerical 

integration  of  the  rotating-frame,  coupled  equations  of  motion: 

^2L  =  _M>l_ym  AB(t) 
dt  T2  y 

dMv  , ,  My 

— T- =  yM  x  AB(t)  -  —  -  yMzB, 
dt  1 2 

dM^  =  yM  B,-Mz+Mo 
dt  1  y  1  T,  T, 

AB(t)  -  Bres  -  (B0  +  Bm  cos(comt)) 

in  which  B0  is  the  external  field,  Bi=y©i  is  the  microwave  field  (y,  the  magnetogyric  ratio  for  the 
electron  spin),  Bres  is  center  of  the  resonance  line  at  the  operative  microwave  frequency,  Bm  is 


the  modulation  amplitude,  com  is  the  modulation  angular  frequency,  and  T i  and  T2  are 
respectively  the  phenomenological  longitudinal  and  transverse  relaxation  times.  The  equations 
of  motion  are  integrated  for  a  series  of  resonance  offsets  (B0-Bres),  and  the  inhomogeneous 
resonance  by  weighting  these  responses  according  to  the  intrinsic  overall  line  shape.  Through 
this  integration  we  obtain  the  time-varying  response  of  the  absorption  (My)  and  dispersion  (Mx) 
signals  at  each  sampled  B0.  By  transformation  of  these  time-domain  signals  the  signal 
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Figure  8.  94.5  GHz  EPR  spectra  of  BDPA  in  polystyrene  (ambient  temperature),  as  function  of 
microwave  power  and  field  modulation  frequency  (amplitude  1G).  Absorption  and  Dispersion  signals  are 
shown;  the  latter  are  registered  in  phase  and  in  quadrature  with  the  field  modulation.  Dispersion  phase  is 
inverted  relative  to  figure  9. 

components  at  harmonics  of  com  are  determined  that  are  in  phase  and  in  phase-quadrature  with 
the  field  modulation.  This  approach  is  not  computationally  efficient,  but  is  very  transparent  in  its 
mathematics  and  performance.  Computations  are  carried  out  on  SGI  Octane  or  Origin  200/2000 
workstations.  (Further  detail  and  C++  code  are  provided  in  Appendix  IV). 

The  primary  results  of  this  work  include: 

•  An  experimental  survey  of  power  and  modulation  frequency  effects  on  the  BDPA/PS 
sample;  exemplary  results  are  shown  in  Figure  8.  For  clarity,  only  the  1st  harmonic  in-phase 
absorption  and  in-phase/phase-quadrature  dispersion  signals  are  shown.  At  the  lower  power 
illustrated,  the  spectral  properties  are  independent  of  modulation  frequency  and  correspond 
essentially  to  the  linear,  slow-passage  regime.  At  the  higher,  saturating  power  level,  the 
absorption  signal  remains  at  modest  amplitude,  but  both  dispersion  signals  grow  dramatically. 
The  phase-quadrature  component  has  the  shape  of  an  undifferentiated  absorption  line,  retains  an 
“authentic”  line  width,  and  exhibits  a  amplitude  maximum  when  com«l/Ti.  The  in-phase  signal 
has  a  dispersion  derivative  line  shape  at  low  com  that  evolves  to  a  power-broadened 
undifferentiated  absorption  shape  with  increasing  com.  The  crossover  point  (zero  amplitude  at  the 
line  center)  occurs  also  when  com  is  close  to  l/Tj. 
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Figure  9.  EPR  signals  simulated  by  numerical  integration  of  the  Bloch  equations.  Dispersion  mode  in  phase 
(solid  line)  and  in'phase-quadrature  (dashed  line)  with  the  magnetic  field  modulation.  Simulations  were  done 
with  the  following  parameters:  Gaussian  linewidth  =  9  G,  T j,  =  18  ps,  Tj  =  lps  ,  Bt  =  0.065  G,  Bm  =  0.5  G. 

•  Simulation  of  the  behavior  of  these  signals  (Figure  9;  only  the  1st  harmonic  dispersion 
signals  are  shown)  using  the  spectral  and  relaxation  parameters  reported  by  Farrar  et  al.[32]. 
These  results  immediately  underscore  the  idea  that  dispersion  rapid  passage  signals  can  provide 
substantially  enhanced  signal  amplitudes.  They  suggest,  moreover,  that  saturating  powers  can  be 
reconciled  with  good  spectral  resolution,  particularly  in  the  case  of  the  quadrature  dispersion 
signal.  This  display  also  has  a  line  shape  that,  as  compared  to  the  in-phase  dispersion  signal,  is 
insensitive  to  com  (and  to  relaxation  properties).  It  therefore  appears  well-suited  to  the 
determination  of  spectroscopic  parameters;  knowledge  of  relaxation  properties  would  be  needed 
to  quantitatively  simulate  the  amplitude,  but  not  the  shape.  This  notion  requires  further  testing  in 
systems,  such  as  the  silica  test  sample,  that  exhibit  more  complicated  line  shapes.  Moreover,  it 
is  important  to  test  the  robustness  of  the  behaviors  described  above  in  particular  with  respect  to 
variations  in  relaxation  properties.  This  testing  will  be  done  through  extensive  numerical 
simulations. 

•  An  experimental  study  of  the  fused  silica  E’  center  as  a  function  of  power  and 
modulation  frequency.  The  sample  used  for  these  experiments  was  a  tube  of  dimensions  0.87 
mm  OD,  0.7  mm  ID,  subjected  to  6xl05  rad  60Co  y-irradiation.  Selected  results  are  shown  in  the 
figures  on  the  next  inserted  pages.  Each  page  exhibits  a  series  of  spectra  obtained  at  the 
indicated  modulation  frequency  and  modulation  amplitude  of  3  G.  Each  spectrum  entails  a  sweep 
of  200  G  at  a  rate  of  120  G/min  and  a  time-cinstant  of  0.3  sec.  The  vertical  scales  on  each  page 
indicate  the  relative  amplitudes  of  the  signals  on  all  of  the  pages.  Three  spectral  displays  are 
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shown,  namely  the  (low-power)  absorption  and  dispersion  signals,  and  the  dispersion  component 
detected  n/2  out  of  phase  with  the  driving  field  modulation.  The  data  on  each  page  were 
obtained  at  a  particular  power  level,  as  indicated.  The  fundamental  result  that  emerges  from  this 
study  is  that  the  behavior  observed  is  analogous  to  that  of  the  BDPA/PS  system  provided  only 
that  the  weighting  of  spin  packets  is  adjusted  to  reflect  the  different  inhomogeneous  line  shapes  - 
a  simple  Gaussian  for  BDPA/PS,  a  Zeeman  powder  pattern  for  the  E’  center  in  quartz. 
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Appendix  I 


Schematics  and  Specification  of  the  94.9  GHz  Microwave  Bridge 


MICROWAVE  HETERODYNE  TWO-CHANNEL  BRIDGE 
OF  THE  COHERENT  W-BAND  EPR  SPECTROMETER 
FOR  PULSED  AND  CW  EXPERIMENTS 


Applications 

The  high-frequency  bridge  represents  a  multifunctional  W-band  transmitter-receiver 
module  (TRM)  intended  for  an  EPR  spectrometer  operation.  The  TRM  has  two  coherent 
transmitter  channels  to  form  a  probing  signal  at  the  central  frequency  of  94.9  GHz  with  the 
maximum  total  power  of  280  mW  as  well  as  the  high-sensitive  receiver  with  double  frequency 

conversion  for  an  EPR  signal  detection. 

Being  capable  of  either  pulsed  or  CW  coherent  operation,  the  system  can  be  operated  in 
several  configurations.  If  integrated  into  measuring  equipment  the  TRM  can  be  employed  in 
scientific  research. 

System  Description  and  Functional  Characteristics 

A  simplified  functional  block  diagram  of  the  TRM  is  shown  in  figure  1.  The  system 
consists  of  several  main  blocks:  the  double-channel  transmitter,  the  receiver  with  double 
frequency  conversion,  the  master  oscillator  (STAMO),  and  the  transmission-reception 

circulator.  The  system  performs  in  the  following  manner.  The  STAMO  supplies  coherent 

centimeter-wave  driving  signals.  These  signals  are  transferred  to  W-band  by  an  IMPATT 
Frequency  Multiplier  with  multiplication  factor  as  high  as  13.  The  transfer  is  realized  in  the 
transmitter  channels  and  receiver  heterodyne  1  without  significant  phase  noise  degradation. 

The  transmitter  signal  AM  and  two-level  phase-shift  keying  (2PSK)  are  accomplished  by 
means  of  external  TTL  control  pulses  supplied  to  the  transmitter  channels.  The  power  gam  is 
achieved  due  to  3-stage  CW  IMPATT  Amplifiers.  The  TRM  cdn  be  operated  in  an  External 
Synthesizer  configuration  to  attain  the  best  frequency  stability  and  possibility  of  frequency 
adjustment  Coherent  summation  of  the  outputs  of  two  transmitter  channels  is  carried  out  by 

means  of  an  E-H  Hybrid  Tee. 

The  receiver  is  designed  on  the  basis  of  a  double  frequency  conversion  circuit  and  has 
two  outputs  which  are  90°  out  of  phase.  The  Balanced  Mixer  input  of  the  receiver  is  protectd 
with  a  High-Speed  PIN  Switch  which  provides  maximum  attenuation  while  the  transmitter 
signals  pass.  For  the  transmitter  noise  to  be  suppressed  during  a  pause  between  the  transmitter 
pulses  another  High-Speed  PIN  Switch  is  included  at  the  transmitter  output. 
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STAMO 

As  mentioned  above,  the  STAMO  provides  the  coherent  centimeter-wave  driving  signals 
for  both  the  transmitter  and  receiver  and  can  perform  in  two  configurations. 

1.  Internal  Stabilized  Oscillator  Configuration.  In  this  configuration,  the  TRM 
performs  at  a  constant  frequency  of  94.9  GHz.  The  system  output  frequency  stability  is 
determined  by  the  stability  of  a  Transistor  Oscillator  at  7.16  GHz  stabilized  with  a  high-Q 
Dielectric  Resonator. 

2.  External  Synthesizer  Configuration.  In  this  configuration,  the  reference  frequency 
signal  (7.16±0.004)  GHz  is  sampled  from  an  External  Synthesizer.  The  operating  output 
frequency  in  the  range  of  (94.9±0.05)  GHz  and  output  frequency  stability  are  determined  by  the 

External  Synthesizer  signal  parameters. 

As  seen  from  figure  1,  the  STAMO  forms  the  following  signals: 

1)  the  signal  for  driving  the  transmitter  channel  1. 

The  central  frequency  is  7.3  GHz. 

2)  the  signal  for  driving  the  transmitter  channel  2. 

The  central  frequency  is  7.3  GHz.  There  is  placed  a  manually  adjustable  phase  shifter  in 
the  channel. 

3)  the  signal  for  driving  the  receiver  heterodyne  1 . 

The  central  frequency  is  7.16  GHz.  There  is  placed  a  manually  adjustable  phase  shifter  in 
the  channel. 

4)  the  signal  for  driving  the  receiver  heterodyne  2. 

The  central  frequency  is  1.82  GHz. 

5)  the  signal  for  testing  the  transmitter  reference  frequency. 

An  interconnection  scheme  of  the  5  f  AMO  is  shown  in  figure  2. 

Transmitter 

The  transmitter  comprises  two  identical  coherent  channels.  Each  transmitter  channel  is  a 
converter-amplifier  appliance  which  enables  transferring  the  STAMO  signal  at  7.3  GHz  into 
the  signal  at  94.9  GHz  without  significant  noise  phase  degradation  and  amplifying  the  last  one 
up  to  ~  300  mW.  Efficient  frequency  conversion  is  parried  out  byx  an  IMP  ATT  CW  Frequency 
Multiplier  with  a  multiplication  factor  as  high  as  13  and  conversion  losses  no  more  than  16  dB. 
The  transmitter  channel  power  gain  is  carried  out  by  means  of  a  3-stage  IMP  ATT  CW  Amplifier 
with  17  dB  gain. 

For  output  power  control,  included  at  the  output  of  each  transmitter  channel  is  a  current- 
controlled  PIN  Attenuator  with  power  suppression  of  40  dB.  The  output  signal  AM  for  each 
transmitter  channel  is  realized  by  2-stage  High-Speed  PIN  Switches  1,2  with  isolation  of  70  dB 
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and  switching  time  about  of  10  ns.  Such  a  value  of  isolation  is  attained  due  to  the  architecture  of 
these  2-stage  Switches  both  consisting  of  two  High-Speed  PIN  Switches  connected  in  series. 

The  output  signal  [0-rc]  2PSK  for  each  transmitter  channel  is  carried  out  by  means  of  a 
High-Speed  PIN  Two  Level  Phase  Shifter.  Coherent  summation  of  the  outputs  of  two  transmitter 
channels  is  provided  by  the  E-H  Hybrid  Tee,  as  mentioned  above. 

To  suppress  the  noise  of  amplifying  stages  during  a  pause  between  the  pulses  High-Speed 
PIN  Switch  3  with  suppression  no  less  than  30  dB  is  included  at  the  transmitter  output. 

At  the  transmission-reception  Circulator  output  each  of  th^  transmitter  channels  supplies 
an  output  power  no  less  than  70  mW,  the  maximum  total  power  cf  both  channels  being  about  of 
280  mW  at  the  central  frequency. 

Receiver 

The  receiver,  is  designed  on  the  basis  of  a  double  frequency  conversion  circuit.  The 
STAMO  sends  two  signals  to  the  receiver:  the  driving  signal  for  the  receiver  heterodyne  1  at 
7.16  GHz  and  that  for  the  receiver  heterodyne  2  at  1.82  GHz.  An  IMPATT  CW  Frequency 
Multiplier  converts  the  7,16  GHz  signal  into  the  93.08  GHz  signal  with  a  multiplication  factor 
as  high  as  13  in  heterodyne  1  track.  For  the  Multiplier  AM  noise  to  be  suppressed  a  Band-Pass 
Filter  is  placed  within  heterodyne  1  track. 

The  heterodyne  power  enters  the  Balanced  Mixer  input  through  High-Speed  PIN  Switch  4 
due  to  which  it  is  modulated.  The  heterodyne  power  modulation  is  injection-locked  with  the 
transmitter  outputs  and  thus  the  linear  behaviour  of  the  receiver  is  provided. 

An  EPR  signal  enters  the  Balanced  Mixer  signal  input  through  the  transmission-reception 
Circulator.  The  signal  input  of  the  Balanced  Mixer  is  protected  with  High-Speed  PIN  Switch  5. 
A  series  of  the  transmitter  signals  passing,  the  Switch  provides  isolation  of  40  dB. 

The  IF  signal  at  1.82  GHz  is  amplified  by  a  Low-Noise  Preamplifier,  divided  by  two 
equal  parts  and  they  enter  the  inputs  of  Mixers.  The  IF  signal  bandwidth  is  1 50  MHz  and  is 
determined  by  a  bandwidth  of  the  Band-Pass  Filter.  After  the  second  frequency  conversion  the 
receiver  has  two  outputs  which  are  out  of  phase  at  the  zero  LF.  The  receiver  provides  a  noise 
figure  no  more  than  12.5  dB  and  power  gain  20  dB  with  a  through  frequency  band  no  less  than 
150  MHz. 

An  interconnection  scheme  of  the  TRM  is  shown  in  figure  3. 
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Product  Design 

The  product  is  fabricated  as  a  stand-alone  unit  with  dimensions  470x300x  120  cm3.  Led 
out  on  the  back  panel  is  the  unit  waveguide  output  connected  to  the  EPR  spectrometer 
waveguide.  The  waveguide  size  is  WR-1 0.  The  waveguide  flangevis  UG-387/U. 

On  the  front  panel,  as  seen  from  figure  4,  there  are  allocated  all  of  the  controls  of  the 
unit  as  well  as  the  supply  voltage  connector,  the  latter  being  arranged  in  the  left  lower  comer 
of  the  panel.  The  microwave  coaxial  connectors  OSM(SMA)  Jack  arc  utilized  for  control  signals 
to  be  supplied. 

EPR  SIGNAL  OUTPUT  1  -  second  IF  output  of  channel  1  of  the  receiver. 

EPR  SIGNAL  OUTPUT  2  -  second  IF  output  of  channel  2  of  the  receiver  which  has  90° 

I 

of  phase  shift  with  respect  to  OUTPUT  1. 

R  PS  ADJ  -  manual  phase  adjustment  for  the  receiver  heterodyne  1. 

T  PS  ADJ  -  manual  phase  adjustment  for  the  transmitter  channel  2. 

RJF  TEST  -  (the  left  side  of  the  panel)  the  transmitter  reference  frequency  testing  at 

7.3  GHz. 

EXT-INT  -  operating  configuration  switch.  "EXT"  position  corresponds  to  the  External 
Synthesizer  configuration  whereas  "1NT"  position  Corresponds  to  the  Internal  Stabilized 
Oscillator  configuration. 

RF  EXT  IN  -  External  Frequency  Synthesizer  input  at  (7.16+0.004)  GHz. 

TST  ON  -  thermostat  system  switch  for  the  Transistor  DRO  at  7. 1 6  GHz. 

REC  CONT  TTL  G ATE  -  the  receiver  control  connector. 

While  applying  TTL  pulses,  PIN  Switch  5  included  at  the  receiver  input  is  on  provided 
these  pulses  are  not  injection-locked  to  the  pulses  of  the  transmitter  channels. 

TRANS  CONT  -  series  of  the  transmitter  control  connectors: 

AM  CHAN-1  TTL  TRIG  -  connector  for  the  transmitter  AM  to  be  triggered  by  TTL 
control  pulses  supplied  to  channel  1. 

AM  CHAN-2  TTL  TRIG  -  connector  for  the  transmitter  AM  to  be  triggered  by  TTL 
control  pulses  supplied  to  channel  2. 

2PSK  CHAN-1  TTL  TRIG  -  connector  for  the  transmitter  [0-ji]  two-level  phase-shift 
keying  (2PSK)  to  be  triggered  by  TTL  control  pulses  supplied  to  channel  1 . 

2PSK  CHAN-2  TTL  TRIG  -  connector  for  the  transmitter  [0-rr]  2PSK.  to  be  triggered  by 

TTL  control  pulses  supplied  to  channel  2. 

ATT  CONTROL  DC  CHAN-1  -  connector  for  the  transmitter  output  to  be  attenuated  by 
control  DC  supplied  to  channel  1 . 
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A1  .  CONTROL  DC  CHAN-2  -  connector  for  the  transmitter  output  to  be  attenuated  by 
control  DC  supplied  to  channel  2. 

OPER  MODE 

PULSE-CW  -  the  TRM  operating  mode  switch.  "PULSE"  position  corresponds  to  the 
pulsed  mode  whereas  "CW"  position  corresponds  to  the  CW  mode  of  operation. 

MEASUR-TUNING  -  operating  submode  switch  which  allows  the  TRM  to  pass  from 

tuning  condition  to  measuring  one. 

Pulsed  mode 

In  the  pulsed  mode,  the  TRM  supplies  AM  modulated  microwave  output.  For  each 
transmitter  channel  the  microwave  pulse  length  and  period  of  pulse  repetition  are  determined  by 
the  parameters  of  the  TTL  control  pulses  supplied  to  the  input  connectors  TTL  TRIG  AM 
CHAN- 1  and  TTL  TRIG  AM  CHAN-2. 

Turning  on  and  off  the  receiver  is  attained  by  the  TTL  control  pulses  supplied  to  the 
connector  TTL  GATE.  The  transmitter  signal  AM  and  the  receiver  switching  are  carried  out  by 
means  of  five  High-Speed  PIN  Switches  (see  block  diagram,  figure  1). 

High-Speed  PIN  Switch  1  provides  the  transmitter  channel  1  output  AM  with  isolation  of 

70  dB. 

High-Speed  PIN  Switch  2  provides  the  transmitter  channel  2  output  AM  with  isolation  of 

70  dB. 

High-Speed  PIN  Switch  3  with  isolation  of  30  dB  introduces  an  additional  depth  of 
modulation  of  the  transmitter  output  and  suppresses  the  CW  IMPATT  Amplifier  noise  of  the 
transmitter  channels  1  and  2  during  a  pause  between  microwave  power  pulses  as  well. 

The  transmitter  power  flowing,  High-Speed  PIN  Switch  4  turns  off  the  receiver 
heterodyne  1  power,  which  permits  the  receiver  to  perform  in  linear  conditions. 

High-Speed  PIN  Switch  5  with  isolation  of  40  dB  protects  the  receiver  Balanced  Mixer 

from  the  transmitter  microwave  power  pulses. 

For  the  proper  order  of  switching  on  and  off  the  High-Speed  PIN  Switches  there  is  a 

logic  and  five  drivers  in  the  system.  Figure  5  shows  a  switching  scheme  of  the  High-Speed  PIN 
Switches  as  a  function  of  the  cpntrol  signals  derived  from  the  External  TTL  Pulse  Oscillator  and 

supplied  to  the  control  connectors  TTL  TRIG  AM  and  TTL  GATE. 

"On"  Switch  state  corresponds  to  the  minimum  loss  performance,  while  "Off"  state 

corresponds  to  the  maximum  isolation  performance. 
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The  TTL  high  ("1")  signal  supply  to  the  input  connector  TTL  TRIG  AM  results  in  the 
transmitter  channels  power  being  emitted  while  the  TTL  low  ("0")  signal  supply  causes  their 
power  being  suppressed. . 

The  receiver  is  driven  by  the  control  signal  entering  the  input  connector  TTL  GATE, 
the  TTL  high  ("1")  signal  turning  the  receiver  on  and  TTL  low  ("0")  signal  turning  the  receiver 
off.  However  (as  seen  from  figure  5,  versions  a,  b,  c),  any  transmitter  channel  emitting 
microwave  power  pulses,  the  receiver  is  always  off  irrespective  of  a  signal  at  the  control 
connector  TTI  GATE.  The  receiver  is  on  only  for  a  case  of  the  transmitter  power  being 
suppressed  when  the  TTL  high  signal  enters  the  connector  TTL  GATE  (version  e).  This  prevents 
the  possibility  of  the  Balanced  Mixer  failure. 

While  operating  in  the  pulsed  mode,  there  is  the  possibility  of  2PSK.  the  output  signal  for 
each  transmitter  channel.  The  TTL  low  or  TTL  high  signal  supply  to  the  input  connector  2PSK 
CHAN-1  TTL  TRIG  results  in  180°  of  phase  change  of  the  transmitter  channel  1  output.  The 
TTL  low  or  TTL  high  signal  supply  to  the  input  connector  2PSK  CHAN-2  TTL  TRIG  results  in 
1 80°  of  phase  change  of  the  transmitter  channel  2  output. 

The  output  pulse  power  of  each  transmitter  channel  can  be  attenuated  by  a  control  current 
supply  to  the  input  connectors  ATT  CONTROL  DC  CHAN-1  and  ATT  CONTROL  DC  CHAN-2. 
Figure  6  shows  the  attenuation  as  a  function  or  a  control  current  for  channel  1  and  channel  2. 

The  pulsed  mode  has  two  submodes:  the  measuring  submode  and  tuning  one. 

In  the  measuring  submode,  the  transmitter  emits  micrdwave  power  pulses  which 
duration  and  period  of  repetition  are  determined  by  the  control  pulse  parameters.  After  double 
frequency  conversion  and  addition-subtraction  processing  the  EPR  signal  detected  by  the 
receiver  is  sent  to  the  connectors  EPR  SIGNAL  OUTPUT  1  and  EPR  SIGNAL  OUTPUT  2. 
The  output  signals  at  these  connectors  have  90°  of  phase  shift  with  respect  to  each  other.  In 
the  measuring  submode,  the  transmitter  microwave  pulses  don't  penetrate  into  the  receiver  due  to 

high  isolation  between  the  receiver  and  transmitter. 

Tuning  submode  is  employed  for  the  EPR  spectrometer  resonator  being  tuned.  In  this 
submode,  isolation  between  the  transmitter  and  receiver  is  intentionally  reduced,  which  is 
achieved’ by  switching  High-Speed  PIN  Switch  4  to  "on"  state,  and  as  a  result  the  receiver 
heterodyne  1  power  enters  the  Balanced  Mixer  even  though  for  a  case  of  the  transmitter  pulses 

being  emitted. 

CWMode 

In  the  CW  mode,  the  TRM  supplies  an  output  power  of  5  mW.  There  is  utilized  only 
channel  2  of  the  transmitter  in  this  mode.  Channel  1  of  the  transmitter  is  completely  off. 
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So  as  to  avoid  the  Balanced  Mixer  overload  and  failure  the  current-  controlled  PIN 
Attenuator  included  at  the  output  of  channel  2  introduces  continued  attenuation  of  8  dB.  The 
output  power  control  can  be  accomplised  only  towards  power  level  reduction. 

In  the  CW  mode,  High  Speed  PIN  Switches  2, 3, 4,5  are  in  "on”  state,  which  corresponds 
to  the  minimum  insertion  loss  performance.  The  High-Speed  PIN  Two  Level  Phase  Shifter  of 
channel  2  is  off. 

Operating  Order 

1.  Attach  the  TRM  output  waveguide  flange  to  the  waveguide  track  of  the  EPR 
spectrometer. 

2.  Supply  the  TRM  with  voltages  by  means  of  a  feeding  cable  in  accordance  with  a 
table  of  supply  voltages. 

Attention.  The  voltage  of  +1 2  V  is  allowed  to  be  applied  to  the  supply  voltage 
connector  only  after  supplying  the  voltage  of  minus  12  V.  Turning  the  supply  voltages  off 
should  be  performed  in  a  reverse  order. 

3.  While  operating  in  the  Internal  Stabilized  Oscillator  configuration,  put  the  switch 

1NT-EXT  (the  right  side  of  the  panel)  into  "INT"  position. 

Testing  of  the  transmitter  reference  frequency  of  7300  MIIz  can  be  carried  out  by 
connecting  a  frequency  meter  to  the  connector  RF  TEST.  The  testing  signal  power  at  the 
connector  RF  TEST  is  about  1  mW. 

4.  While  operating  in  the  External  Synthesizer  configuration,  supply  the  External 
Synthesizer  signal  to  the  connector  RF  EXT  IN.  The  signal  frequency  is  within  the  range  (7156- 
7164)  MHz.  The  Synthesizer  power  output  is  of  (5-10)  mW. 

Put  the  switch  INT-EXT  (the  ripht  side  of  the  panel)  into  "EXT"  position. 

5.  While  operating  in*the  pulsed  mode,  put  the  switch  PULSE-CW  into  PULSE 

position. 

While  tuning  the  EPR  spectrometer  resonator,  put  the  switch  MEASUR-TUN1NG  into 
"TUNING"  position. 

While  measuring,  put  the  switch  into  "MEASUR"  position. 

The  control  signals  should  be  supplied  to  the  connectors  REC  CONT  and  TRANS 

CONT  in  accordance  with  a  section:  Pulsed  Mode. 

The  transmitter  channel  2  phase  adjustment  should  be  carried  out  by  turning  the  knob 

T  PS  ADJ. 

The  receiver  heterodyne  l  phase  adjustment  should  be  carried  out  by  turning  the  knob 
R  PS  ADJ. 
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After  double  frequency  conversion  the  EPR  signal  delected  by  the  receiver  enters  the 
connectors  EPR  SIGNAL  OUTPUT  1  and  EPR  SIGNAL  OUTPUT  2.  The  output  signals  at 
these  connectors  have  90°  of  phase  shift  with  respect  to  each  other. 

6.  While  operating  in  the  CW  mode,  put  the  switch  PULSE-CW  into  "CW"  position 
and  the  switch  MEASUR-TUNING  put  into  "MEASUR"  position. 

The  transmitter  channel  2  phase  adjustment  should  be  carried  out ;  by  turning  the  knob 
T  PS  ADJ. 

The  receiver  heterodyne  1  phase  adjustment  should  be  carried  out  by  turning  the  knob 
R  PS  ADJ. 

After  double  frequency  conversion  the  EPR  signal  detected  by  the  receiver  enters  the 
connectors  EPR  SIGNAL  OUTPUT  1  and  EPR  SIGNAL  OUTPUT  2.  The  outputs  at  these 
connectors  have  90°  of  phase  shift  with  respect  to  each  other. 

Duty  conditions 

1.  The  product  is  to  be  incorporated  in  the  EPR  spectrometer  and  employed  for  in-door 
operation  under  standard  climatic  conditions  free  of  mechanical  perturbations: 

•  Ambient  Temperature,  °C . 25+10 

•  Relative  Humidity,  % . . . 70 

2.  The  product  should  be  installed  so  that  a  value  of  an  external  magnetic  field  and  its 
deviation  for  a  single  EPR  experiment  could  not  exceed  1  tnT  and  0.02  mT,  respectively. 

3.  A  1 1  e  n  t  i  o  n.  While  transporting,  storing,  assembling,  and  operating,  it  is 
necessary  to  preclude  entering  any  particles  (especially  magnetic  ones)  into  the  waveguide 
channel  so  as  to  avoid  the  product  failure. 

4.  Warming-up  time  after  switching  on  is  1  hour  (max). 

5.  Continuous  duty  time  is  8  hours. 

6.  Guarantee  period  is  2  years. 
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SPECIFICATIONS 


1.  TRANSMITTER 

PULSE  MODE 

Internal^  Slabilized_  Oscillator  Con[h 
Characteristic,  units  [ 


1 .  Output  frequency  fixed,  GHz  _ 


2.  Output  frequency  instability,  kHz/h 


3.  Output  spectral  linewidth  at  the  level 
-3dB,  kHz _ 


4.  AM  noise  with  offset  from  carrier, 
dB/Hz 
1  kHz 
10  kHz 

100  kHz _ 


5.  Output  power  of  channel  l,  mW _ 


6.  Output  power  of  channel  2,  mW _ 


7.  Output  power  attenuation  range,  dB 
Channel  1 

Channel  2  _ _ 


8.  Maximum  output  power  of  two 
channels,  mW _ _ _ 


9.  Microwave  pulse  duration  at  the  level 
-3dB,  ns 


10.  Pause  duration  between  microwave 
pulses,  ns 

min 

max _ _ _ 

1 1.  Microwave  power  switching  time  at 

the  level  0. 1/0.9,  ns _ . 


12.  Microwave  power  suppression  during 
a  pause  between  pulses,  dB _ _ 


13.  Switching  time  up  to  maximum 

suppression  level,  ns _ _ 

14.  Accuracy  of  180  deg  phase  shift 

keeping  of  2PSK  for  94.9  GHz 
frequency,  deg _ 


15.  Switching  time  for  2PSK,  ns _ 

16.  Output  parasite  AM  level  for  2PSK, 
dB 


uralion _ 

Standards  of 
Performance 
Specifications 


94.90 


150 


continuous 


20 

continuous 


Measured  Value 


94.900 


130 


20 

continuous 


20 

continue'. 
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2.  TRANSMITTER 
PULSE  MODE 


External  Synthesizer  Configuration 


Characteristic,  units  _ 


1 .  Output  frequency,  GHz  _ 


2.  Output  power  of  channel  1,  mW 


3.  Output  power  of  channel  2,  mW 


4.  Maximum  output  power  of  two 
channels,  mW _ 


5.  Microwave  power  suppression 
during  a  pause  between  pulses,  dB 


6.  Accuracy  of  180  deg  phase  shift 
keeping  for  2PSK,  deg _ 


7.  Output  parasite  AM  level  for 
2PSK,  dB 


94.770 


Measured  Value 


94.835  94.900  94.965 


95.030 


3.  TRANSMITTER 


4.  RECEIVER 

Internal  Stabilized  Oscillator  or  External  Synthesizer  C  onfigurations 


Characteristic,  units _ 

1.  Input  signal  frequency,  GHz 


2.  Noise  figure,  dB _ __ 


3.  Power  gain,  EPR  signal  1 
dB 


4.  Through  band,  MHz 


Measured  Value 


94.770 


94.835 


94.900 

11.6 


94.965 

11.8 


95.030 
1 1.7 


5.  TRANSMITTER  -  RECEIVER 

1.  Total  phase  drift  *  5deg/15  minutes. 

2 .  Manual  signal  phase  adjustment 

for  channel  2  .  more  than  360  deg. 

3 .  Manual  signal  phase  adjustment 

for  receiver  heterodyne  1  ■ more  ^an  ^60  deg. 
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6.  SUPPLY  VOLTAGE  AND  DRAIN  CURRENT  PARAMETERS 


Supply  Voltage 

Drain  Current 

Pulsed  Mode 

CW  Mode 

U1  =  +27V±  IV 

11  =  1.82-1.53  A 

11  =0.9  A 

U2  =  +  12V±0.5V 

12  =  2.5  A 

12=  1.97  A 

U3=  -12V  ±  0.5V 

13  =  0.45  A 

13  =  0.42  A 

7.  THE  TTL  CONTROL  PULSE  PARAMETERS 
TTL  high  ("I")  (+2.4  - +4)  V 

TTL  low  ("0")  (0  -  +0.4)  V 

Switching  time  0. 1/0.9  x  V  no  more  than  1  ns 
Load  resistance  50  ohm 


8.  REFERENCE  DATA 

1.  The  output  power  at  the  connector  RF  TEST  is  about  1.6  mW. 

2.  The  input  power  supplied  by  the  External  Synthesizer  to  the  connector 
RF  EXT  IN  is  (5-10)  mW. 

3.  The  output  signal  frequency  of  the  transmitter  is  related  to  the  input  signal 
frequency  of  the  External  Synthesizer  by  the  following  relationship: 

Fout  t(MHz)  =  Finput  syn(MHz)  +  140]  x  13. 
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X2  'RF  EXT  IN' 


Interconnection  schene  of  STAMD  (nodule  B). 


Component 

Designation 

Component  Name 

G1 

Transistor  DRO 

G2 

Reference  Crystal  Oscillator 

ARA2 

Transistor  Power  Amplifier 

W1-5-W2 

IF  Band  Pass  Filter 

01 

Microstrip  Directional  Couple. 

02 

Microstrip  Power  Divider 

FI,  F2 

Manual  Adjust  Coaxial  Phase  Shifter 

El 

SSB  Upconverter 

E2  | 

Double-Output  Mixer 

WV1 

Frequency  Multiplier  x  13 

VI 

Thermostat  Control  Circuit 

V2,  V3 

Regulation  Circuit 

B1*B2 

Switch 

K1 

Coaxial  Switch 

XI 

Voltage  Supply  Connector 

X2+X9 

Connector  OSM  (SMA) 

List  of  Component  to  Interconnection  Scheme  of 
the  STAMO 
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(Z-IX)A^*  (CA)tr 


Interconnection'  schene  of  Tronsnitter-Receiver  (nodule  A). 


Component 

Designation 

Component  Name 

WS1-5-WS19 

Waveguide  Circulator 

WV7-rWV9 

IMP  ATT  CW  Frequency  Multiplier 

WVl-rWV6 

IMPATT  CW  Amplifier 

WURWU7 

I-Iiph-SDeed  PIN  Switch 

High-Speed  PIN  Two  Level  Phase  Shifter 

U1.U2 

PIN  Attenuator 

S1-J-S9 

Driver 

Y1-5-Y9 

Current  Regulator 

W1 

IF  Band  Pass  Filter 

W2-rW5 

Waveguide  Band  Pass  Filter 

El 

Balanced  Mixer 

Q1 

E-H  Hybrid  Tee 

A1 

Low  Noise  Preamplifier 

A2^-A4 

Transistor  Power  Amplifier 

V14V6 

Voltage  Regulator  and  Control  Circuit 

XI 

Voltage  Supply  Connector 

X2-rX5 

Connector  OSM  (SMA) 

X6-rX  1 0 

Connector  OSM  (SMA) 

Bl,  B2 

Switch 

List  of  Component  to  Interconnection  Scheme  of 
the  Transmitter-Receiver  Module  (fig.3). 
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94.9  GHz  EPR  spectrometer  microwave  bridge 

Block  diagram  showing  the  W-band  spectrometer  bridge.  The  high  stability  (ST AMO)  master  frequency  oscillator 
system  lies  in  ol'-'ck  (A)  of  the  diagram.  The  two-channel,  phase-adjustable,  pulsed  transmitter  assembly  lies  in 
block  (B).  Block  (C)  is  the  quadrature  receiver  unit. 
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Voltage 
supply  “ 
connector 


Front  Panel 
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TTL  INPUT 

LEVEL  CDNNECTUR 
]  I  TTL  TRIG  AM" 
1  '  CHANl 

,  77  TTL  TRIG  AM 

*  1  CHANg _ 


'1' 

T 

— > 

Tor'0J 

! 

•v 

_ 

- > 

TorT 

- > 

TTL  TRIG  AM 
CHANl 

TTL  TRIG  AM 
CHANg 

TTL  GATE 


— >  DRIVER  1 
-*  DRIVER  g 
^"DRIVER  3 
-*  DRIVER  4 
DRIVER  5 

-4  DRIVER  T 
->  DRIVER  g 
->  DRIVER  3 
-*  DRIVER  4 


DRIVE 

SIGNAL 

^7TlAM-l 


Hiqh-Spped 


AM-g  . 
AM- 3 ' 
AM-4  ' 
AM-5. 

AM-1  , 
AM-g. 
AM- 3 , 
AM-4  * 


PIN 

n  : 

Switch 

SI  ATE 

i 

'on' 

?. 

'on' 

3 

'on' 

4 

_5_ 

'off' 

'off'* 

■  1 

"on' 

-  2 

'off' 

„  0 

'on' 

'off' 

>  5 

'off' 

'0' 

'1' 

— 5, 

Tor'01 

_ 

TTL  TRIG  AM 
CHANl 

TTL  TRIG  AM 
CHANg 

TTL  GATE 


DRIVER  1 


DRIVER  4 
DRIVER  5 


I  AM- 1 


AM -4 
AM -5 


1 

'off' 

0 

L. 

'on' 

3 

'on' 

4 

'off' 

>  5 

'off' 

'0' 

— ? 

'0' 

— > 

'0' 

— ^ 

TTL  TRIG  AM 
CHANl 

TTL  TRIG  AM 
CHANg 

TTL  GATE 


4  DRIVER  1  AM~MT 
DRIVER  g  -^L^“g 
->  DR  t  VCR  3  3 

->  DRIVER  4  -MVL  4 
-*  DRIVER  5^-^^  5 


1 

'off' 

g 

'off' 

_3_ 

"off' 

4 

'off' 

•  5 

'off' 

'0' 

— => 

'0' 

— ^ 

T 

— =» 

TTL  TRIG  AM 
CHANl 

TTL  TRIG  AM 
CHANg 

TTL  GATE 


1  ^  DRIVER  1 . 

'  „  DRIVER  g  ^ 

DRIVER  3  ^ 

-  -*-»  DRIVER  A 

_~*_-4DRIVER 

Switching  Scheme. 


'off' 


4 

'on' 

5  1 

"on' 
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Insertion  Loss,  dB 


Microwave  power  attenuation  as  a  function  of  a  current  value.(Channel  I) 


Appendix  II 


Lab  VIEW  Code  for 

EPR  Spectrometer  Control,  Data  Acquisition,  and  Manipulation 


epr  program  version  1.5.5.vi 


Labview  Routine  for  EPR  Control  and  Data  Acquisition:  Labview  Code 


-$> 


epr  program  version  1.5.5.vi 
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A  Search  for  High-Spin  Zeeman  Effects  in  the  94.9  GHz  EPR  spectrum  of  Cr3+  in  a-Al203 


94.9  GHz  EPR  Study  of  Cr3+  in  (X-AJ2O3 
D.  A.  Schwartz,  E.  D.  Walter  and  D.  J.  Singel 


Abstract 

We  present  results  of  a  high-frequency  (94.9  GHz)  continuous  wave  (CW)  room- 
temperature  electron  paramagnetic  resonance  (EPR)  study  of  Cr3+  (S=3/2)  in  a  small  neutron- 
irradiated  and  thermally  annealed  single  crystal  of  (X-AI2O3  (sapphire)  grown  by  the  laser  heated 
pedestal  method.  Through  analysis  of  the  angular  variation  patterns  of  the  three  nominally 
allowed  electronic  spin  transitions  of  Cr3+,  we  have  measured  spin-Hamiltonian  parameters  (gx  = 
1.9812,  gy  =  1.9820,  |D|  =  5.737  GHz)  and  have  considered  the  potential  impact  of  high-spin 
Zeeman  terms  (terms  of  the  form  BS3)  on  the  outcome  of  spin-Hamiltonian  parameter 
assessment.  Our  results  indicate  the  marginal  relevance  of  higher  order  terms  to  the  satisfactory 
analysis  of  the  observed  angular  variation  patterns  and  ultimately,  to  the  assessment  of  the  point 
symmetry  of  Cr3+  in  sapphire. 

Introduction 

Previously  we  reported  the  observation,  via  EPR  spectroscopy  at  94.9  GHz,  of  a  number 
of  paramagnetic  centers  in  a  single-crystal  fiber  of  sapphire  (01-AI2O3)  that  had  been  treated  by 
neutron-irradiation  followed  by  tnermal  annealing.  In  addition  to  EPR  signals  attributable  to 
radiation-induced  defects  (color-centers),  we  also  observed  resonances  from  Cp+  and  Fe3+  ions. 
The  identity  of  the  ions  was  established  by  comparing  values  of  spin-Hamiltonian  (SH) 
parameters  that  we  estimated  from  the  94.9  GHz  EPR  spectra  with  values  previously  reported  in 
the  literature. 

We  noted  previously  that  the  EPR  literature  on  Cr3+  in  sapphire,  while  extensive,  is  not 
completely  consistent.  We  detail  previous  measurements  of  sapphire:Cr3+  spin-Hamiltonian 

parameters  in  Table  1.  Manenkov  and  Prokhorov  2]  originally  reported  the  room  temperature 
SH  parameter  values  g,  =  1.984,  gj_  =  1.9867,  and  |D|  =  5.732  GHz  on  the  basis  of  work  at  8.96 
and  11.97  GHz.  Subsequent  work  at  37.86  GHz  by  Zverev  and  Prokhorov  4]  yielded 
(presumably  also  at  room-temperature)  gj  =  1.982,  g±  =  1.979  and  D  =  -5.732  GHz,  while 
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Geusic  [^],  working  at  9.309  and  11.593  GHz  and  room  temperature,  reported  g(  —  2.003,  gx  ~ 
2.00,  and  D  =  5.786  GHz,  and  Zaripov  and  Shamonin,  working  at  9.58  GHz,  reported  g,  «  gx  = 
1.98  and  |D|  =  5.70  GHz  [6.  7],  also  at  room  temperature.  In  later  work,  several  other  groups  [8> 
9]  determined  SH  parameter  values  that  were  more  like  the  values  Manenkov  and  Prokhorov 
initially  reported  in  1955.  Still  others  D°],  t11!  cited  die  values  reported  by  Manenkov  and 
Prokhorov. 

Klein  and  co-workers  P]  detailed  the  effect  of  temperature  on  the  EPR  parameters  of 
Cr3+  in  ruby.  They  state  that  within  the  accuracy  of  their  measurements  there  was  no 
temperature  dependence  of  the  g-values,  however,  their  Figure  2  seems  to  show  otherwise  and 
the  data  therein  may  support  an  alternate  explanation  of  the  temperature-dependence  of  the  g- 
values.  They  observed  |D|,  however,  to  increase  from  5.726  GHz  to  5.816  GHz  over  the 
temperature  range  88  K  to  575  K.  Thus,  the  variations  in  the  D  values  reported  previously  might 
be  taken  to  indicate  variations  in  the  experimental  temperatures  of  the  above  studies.  These 
results,  however,  suggest  that  variation  in  observed  |D|  values  due  to  possible  temperature 
variation  (-0.24  MHz/K  near  RT)  cannot  completely  explain  the  range  of  previously  reported  D 
values.  The  variety  of  previously  reported  g-values,  too,  cannot  be  explained  on  the  basis  of  this 
variable  temperature  study. 

The  lack  of  consistency  in  these  reported  g-values  may,  of  course,  simply  reflect 
limitations  present  on  the  quantitative  analysis  of  comolsx  EPR  spectra  at  the  time  of  the  early 
studies.  A  detailed  study,  with  extensive  numerical  analysis,  was  reported  in  1978  by  Chang,  et 

a/., [12]  who  found:  g(  =  1.9817  and  gx  =  1.9819,  and  D  =  -5.746  GHz. 

A  consistent  feature  of  all  previous  EPR  studies  of  sapphire:Cr3+  is  the  use  of  a 
“conventional”  spin-Hamiltonian  (CSH)  for  summarizing  the  experimental  results.  This  CSH, 
an  effective  Hamiltonian  operating  on  the  ground  state  crystal  field  manifold,  derives  from  a 
perturbation  treatment  of  the  spin-orbit  coupling  and  electronic  Zeeman  interaction  within  the 

zero-order  crystal  field  states^],  and  contains  only  three  parameters,  g(,  gx,  D.  It  does  not 
contain  the  HSZ  (“high-spin  Zeeman”)  terms  -  terms  linear  in  applied  magnetic  field  and  tri- 
linear  in  spin  operators  -  that  are  included  for  S=3/2  centers  in  the  GSH  (“generalized”  spin 

Hamiltonian)  treatments  presented  by  numerous  authors  (see,  e.g.,  RoitsinD^],  McGavin,  et 
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a/.  [15]  and  references  cited  in  loc.).  While  HSZ  terms  are  often  dismissed  a  priori  as  negligible, 
Geusic,  specifically  called  attention  to  their  possible  importance  in  his  1956  study  of  the  EPR  of 

Cr3+  in  sapphiret^l. 

As  operating  frequencies  are  increased,  the  HSZ  interactions  could  gain  in  significance 
relative  to  the  field-independent  fine  structure  interaction.  The  capacity  to  measure  such 
interactions  could  be  an  important  singular  achievement  of  the  emerging  held  of  “high- 
frequency”  (£100  GHz)  EPR  spectroscopy  (see,  e.g.,  Buckmaster  et  One  benefit  of 

such  measurements  relates  to  the  ability  of  SH  terms  supralinear  in  spin  operators  to  reflect  the 

point  group  of  the  paramagnetic  ion,  as  discussed  in  detail  by  RoitsinP^]  and  by  GrachevP?]. 

In  the  conventional  approach  the  32  crystallographic  point  groups  subdivide  into  three 
categories,  characterized  by  isotropic,  axial,  and  rhombic  tensorial  interactions.  The  applicable 
GSH’s,  however,  subdivide  in  as  many  as  eleven  distinct  symmetry  categories,  or 

“aggregates”^].  Thus,  for  example,  systems  of  C 3,  Civ,  and  C 3h  symmetries  possess  distinct 
spin-Hamiltonians;  if  the  distinguishing  terms  in  the  SH  are  of  sufficient  magnitude,  such 
systems  could  be  distinguished  by  magnetic  resonance  spectroscopy.  Such  discriminating  power 
would  be  of  great  value,  particularly  for  any  variation  from  idealized  uniaxial  symmetry  in  the 
environments  of  paramagnetic  centers,  e.g.,  in  bioinorganic  complexes. 

With  the  goals  of  precisely  characterizing  the  angular  dependence  of  the  EPR  of  Cr3+  in 
sapphire,  and  of  assessing  the  significance  of  the  HSZ  interactions  at  high  field  in  this  system, 
we  have  carried  out  a  study  of  the  orientation  dependence  of  the  room-temperature  EPR  spectra 
of  sapphire:Cr3+  at  94.9  GHz.  We  have  measured  Zeeman  (g,,  gj_)  parameters  and  the  fine 
structure  parameter  (|D|)  from  the  angular  variations  of  resonance  fields  under  the  assumption  of 
uniaxial  symmetry.  We  have  also  analyzed  the  observed  orientation  dependence  of  resonance 
fields  for  possible  HSZ  effects  consistent  with  various  point  symmetries  possible  in  this  system. 
In  light  of  the  present  results,  we  summarize  some  prerequisites  for  successful  observation  of 
HSZ  effects. 

Theory 

The  CSH  for  an  S=3/2  center  located  at  a  site  of  axial  symmetry  in  a  crystal  can  be 
expressed  as: 
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S  =  3/2,  (1) 


<K=  +  Im-bI  lB°l  {gu  cos©  Sz  +  g±  sin©  SJ  +  D  Sz2; 

in  which  pB  is  the  Bohr  magneton,  B0  is  the  applied  magnetic  field,  and  g,,  g±  and  D  are 
parameters  that  gauge  the  Zeeman  (g,,  gj_)  and  fine  structure  (D)  interactions.  In  the  axis  system 
implied  by  equation  (1)  z  is  the  unique  symmetry  axis  in  the  crystal,  and  0  is  the  angle  between 
B0  and  the  z  axis.  For  paramagnetic  centers  with  lower  than  axial  symmetry,  the  Zeeman  portion 
of  the  spin  hamiltonian  is  parameterized  by  gxx,  gyy,  and  gz*,  and  the  fine  structure  interaction  is 
commonly  parameterized  by  D  =  3Dzz/2  and  E  =  (Dxx-Dyy)/2  to  give 


0-C=  +  ||iB|  |B0|  {g^  cos©  Sz  +  gyy  sin©  sin<j)  Sy  +  gxx  sin©  cos<j>  Sx} 

+  DSZ2  +  E  (Sx2  -  Sy2)  (2) 


where  <|>  is  the  angular  displacement  of  B0  from  the  x  axis  about  the  z  axis,  and  the  cartesian  axes 


correspond  to  the  principal  axes  of  the  Zeeman  and  fine  structure  terms. 


Several  authors  have  presented  schemes  for  constructing  more  general  SH’s  (e.g.,  refs  3a, 
4-6  of  Roitsint14!,  refs  1-16  of  Grachev!47],  refs  4  and  5  of  McGavin,  et  a/.H5], 
GeschwindH  8].  This  development  was  initiated  by  Koster  and  Statzl4^],  and  has  been  discussed 

in  detail  by  Roitsinl44],  Buckmaster,  et  20]}  Grachevt^7],  and  McGavin,  et  The 

GSH’s  presented  by  these  authors  are  designed  to  contain  all  terms,  for  a  given  value  of  S,  that 
are  consistent  with  the  symmetry  properties  of  a  crystallographic  point,  and  thus  can  include 
additional  terms  not  present  in  the  CSH  (eq.  (1)). 

For  an  S=3/2  system,  the  additional  terms,  TC  Bs3  --  which  are  linear  in  B0  and  tri-linear  in 

S  [21, 22] ..  are  illustrated  most  lucidly  in  the  formulation  ofRoitsinD4!: 


(3) 
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m 


which:  (Bq  is  a  Cartesian  component  of  the  applied  magnetic  field, 


B™  and  C™  gauge 


the  SH  terms,  and  0“  and  Q™  are  tri-linear  spin  operators  expressed  in  irreducible  tensorial 

form  (tabulated  in  Appendix  V  of  Alt’shuler  and  Kozirevt23]).  As  expressed  in  equation  (3), 
3-fes3  contains  21  adjustable  parameters. 

A  reconstruction  of  the  GSH  has  been  given  by  GrachevD^]  (see  also  Alt’shuler  and 

Kozirev  t23!,  pp.85ff,  and  Grant  and  Strandbergt24!).  This  proceeded  by  coupling  the  field  and 
spin  vector  operators,  decomposing  the  products  according  to  irreducible  tensorial  forms  and 

creating  tesseral  linear  combinations  of  these  tensor  operators  C23>  2^L 

McGavin,  et  al.,  following  this  approach,  give  the  following  terms  for  addition  to  the 
CSH  of  equation  (1): 


O-ChsZ  —  Se  IM-bI  B0 


ZB2>mU2,m 

m=-2 


+  Ib4, 

m=-4 


m^4,m 


(4) 


The  B/>m  parameters  in  Tfosz  are  unitless  variables  indicating  the  size  of  the  interaction  and  the 
U’s  are  operators  that  are  products  of  Hx,  Hy,  Hz,  Sx,  Sy,  and  Sz,  as  detailed  in  McGavin,  et 

a/.  [15]  We  have  dropped  superscripts  for  B’s  since  we  are  only  considering  terms  of  the  form 
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BS3  (i.e.,  our  terms  Bn  m  and  Bn  m  above  are  analogous  to  symbols  B2’ m  and  B4’  m  of 
McGavin;  't  a/.  [15],  eq.  23). 

An  important  advantage  of  this  approach  to  reconstructing  the  GSH  is  that  it  effects  the 
removal  of  certain  redundant  parameters  from  the  SH  of  RoitsinBTi.  The  number  of  terms  BS3 
is  reduced  from  21  in  equation  (3)  to  14  in  equation  (4).  Finally,  as  pointed  out  by  Grachev, 
these  terms  constitute  an  orthogonal  set. 

The  number  of  terms  in  this  qfr,sz  can  be  further  reduced  by  point  symmetry  constraints. 
Each  term  in  the  SH  applicable  to  a  paramagnetic  center  must  transform  as  the  identical  (totally 
symmetric)  representation  of  the  point  group  of  the  center.  To  illustrate  this  adaptation  in  some 
generality  we  consider  C 3  and  all  other  trigonal  groups,  as  well  as  the  hexagonal  groups  with  a 
three-fold  major  proper  rotation  axis.  In  doing  so,  we  aim  to  characterize  what  is  offered  by 
including  these  terms. 
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To  induce  irreducible  representations  of  point  groups  upon  the  SH  terms  in  (4)  we  take 
advantage  of  the  fact  that  they  transform  like  the  tesseral  and  sectoral  harmonic  functions.  We 
therefore  simply  test  for  invariance  of  the  tesseral  and  sectoral  harmonic  functions  with  respect 
to  the  operations  of  a  point  group.  The  test  proceeds  by  considering  the  point  group  symmetry 
operations  St.  These  map  a  general  point  P  in  Cartesian  space,  with  coordinates  (x,y,z),  to  the 
point  P'  with  coordinates  (x',y\z').  Consider  a  function,/ (exemplified  by  a  tesseral  or  sectoral 
harmonic  function)  which  is  expressed  as  /x,y,z).  We  construct  operators-on-functions,  0$, 
corresponding  to  each  of  the  point  group  symmetry  operations  9t  such  that  Og /=/',  where/'  is 
a  new  function  that  assigns  to  a  general  point  P'  the  same  value  as  the  function/assigns  to  P  (see 

Schonlandt27!,  ch.  4).  If  Og/  =  /in  general,  for  all  St  in  a  point  group,  then  /  transforms 
according  to  the  identical  representation  of  that  point  group.  A  spherical  tensor  operator 
analogous  to  the  function  /  (analogous  in  the  sense  of  operator  equivalence  to  the  symmetry 
adapted  Cartesian  function)  then  belongs  to  the  spin-Hamiltonian  appropriate  to  the  point  group 
whose  operations  9t  were  considered. 

The  GSH  for  S=3/2  systems  becomes 

3-C  =  ^  +  ^Hhsz 

where  Q-Q.  is  the  CSH  appearing  equation  (1).  The  high-order  terms  of  equation  (4)  which 
survive  the  Cj  symmetry  operations  are 

^HSZ ;(Cj)  =  +Se  IM-bI  So  (®2,0^2,0  “*"^4, 0^4,0  “**®4,3^4,3  ®4,-3^4,-3  )  (6) 

where  ge  is  the  free  electron  g  value. 

Table  2  shows  the  allowed  terms  from  equation  (4)  for  the  several  point  groups  we  are 
considering.  This  tabulation  illustrates  the  great  advantage  of  this  derivation  of  the  SH:  it  allows 
distinctions  between  point  groups  that  might  otherwise  be  lumped  into  the  “axial”  set  by  analysis 
employing  only  the  CSH  of  equation  (1). 

Materials  and  Instrumentation 

The  single  crystal  of  0C-AI2O3  studied  in  the  present  work  was  fabricated  in  the  laboratory 
of  Professor  William  Yen  (University  of  Georgia)  by  the  laser-heated  pedestal-growth  method 
from  materials  of  99.999%  purityt2^].  This  sample  was  neutron  irradiated,  and  then  thermally 
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annealed  at  -600K  by  Professor  Ilmo  Sildos  (University  of  Tartu,  Estonia).  We  obtained  this 
crystal  -  a  right  circular  cylinder  of  ~0.5  mm  diameter  -  from  Professor  Alexander  Rebane  of  the 
Montana  State  University  Department  of  Physics.  This  crystal  was  found  to  contain  detectable 
amounts  of  Cr3+  and  Fe3+  by  EPR  spectroscopy  at  94.9  GHzt29].  Our  instrumentation,  described 

in  more  detail  elsewheret29!,  includes  a  94.9  GHz  microwave  bridge  and  cavity  obtained  from 
the  Physico-Technical  Institute  of  Donetsk,  Ukraine,  and  a  sweepable  Oxford  Teslatron 
superconducting  magnet.  The  cavity  apparatus  provides  for  frill  rotation  of  the  sample  about  the 
cavity  axis,  and  in  our  instrument  this  axis  may  be  oriented  perpendicular  to  the  applied  static 
magnetic  field  (B0). 

Experiments  and  Results 

CW  EPR  spectra  were  obtained  at  94.9  GHz  as  a  function  of  the  crystal’s  orientation  in 
the  applied  magnetic  field  B0.  Measured  resonance  fields  due  to  Cr3*  are  shown  in  Figure  1  for 
data  obtained  from  a  plane  containing  the  z  axis.  The  sample  was  positioned  for  rotation  by  the 
standard  method:  it  was  placed  in  the  magnetic  field  with  its  cylinder  axis  (nominally  the 
crystallographic  b  axis)  approximately  perpendicular  to  the  applied  magnetic  field  and  nominally 
parallel  to  the  rotation  axis.  The  sample  was  then  oriented  further  by  rotation  about  the  cylinder 
axis  and  another  axis  perpendicular  to  B0  until  the  smallest  value  of  a  resonance  field  Bres  was 

observed  for  any  of  the  strong  high-field  Cr3+  transitions  of  this  sample  (c./  Schwartz,  et  al.  t29]). 
As  a  point  of  departure,  rotation  of  the  sample  proceeded  from  this  point  about  the  cylinder  axis 
with  EPR  spectra  and  resonance  fields  for  the  three  high-field  S=3/2  transitions  being  obtained  at 
intervals  of  A0=O.O4  n  (=  8°)  or  less.  This  rotation  corresponds  (nominally)  to  a  rotation  of  B0  in 
the  crystal’s  xz  plane  where  0  =  0  corresponds  (nominally)  to  B0  ||  c  axis.  Resonance  fields  were 
also  measured  from  other  orientations,  out  of  the  xz  plane  by  adjustment  of  the  microwave 
coupling  adjustment  (which  effects  sample  rotation  about  another  axis  1  B0).  Substantial  mis- 
adjustment  of  microwave  coupling  did  not  interfere  with  the  unambiguous  observation  of  the 
strong  Cr3+  resonances  in  this  sample. 

SH  parameter  optimizations  based  on  observed  AVP’s  (angular  variation  patterns)  of 
resonance  fields  were  conducted  using  the  EPR-NMR  magnetic  resonance  simulation  and 

analysis  software  packaged].  This  software  provides  for  exact  diagonalization  of  a  large 
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variety  of  SH’s  applicable  to  electron  and  nuclear  magnetic  resonances.  Hamiltonians  amenable 
to  analysis  with  this  program  include  general  SH’s  applicable  to  S=3/2  systems  (e.g.  eq.  (5)). 

Analysis  and  Discussion 

Analysis  of  observed  resonance  field  AVP’s  commenced  by  fitting  the  parameters  of  the 
SH  of  eq.  (1)  to  resonance  field  data  in  Figure  1.  The  results  of  this  and  all  other  parameter 
optimizations  were  produced  with  the  computer  program  EPR-NMR  and  are  tabulated  in  Table 
3.  With  linewidths  ranging  from  10-40  Gauss,  the  size  of  the  root  mean  square  deviation 
(RMSD)  between  calculated  and  observed  resonance  fields  for  this  initial  analysis  (24.86  G)  is 
not  unexpected.  Nevertheless,  we  are  interested  in  the  possibility  that  BS3  effects  are  discernible 

in  the  observed  resonance  fields  even  though  these  effects  may  be  very  smallt16!.  Accordingly, 
we  shall  explore  several  explanations  for  the  deviations  between  calculated  and  observed 
resonance  fields  for  this  initial  parameter  determination.  The  deviations,  if  properly  accounted 
for,  may  allow  HSZ  effects  in  sapphire:Cr3+  to  be  quantified.  Contributions  to  the  RMSD 
between  observed  resonance  fields  and  those  calculated  from  the  CSH  include:  1)  effects  on  the 
positions  of  resonance  fields  (including,  but  not  limited  to  effects  on  the  angular  positions  of  the 
extrema  of  the  resonance  field  tracks)  due  to  HSZ  effects;  2)  improper  positioning  of  the  sample 
in  the  magnetic  field  such  that  sample  rotation  was  carried  out  with  some  systematic  error  in  the 
determination  of  the  angular  sample  position;  3)  some  combination  of  1)  and  2);  and  4)  non- 
systematic  errors  in  the  determination  of  sample  angular  positions. 

We  now  consider  the  impact  of  various  Bn>mUn,m  terms  on  the  quality  ox  the  SH  parameter 
fits  to  uncorrected  resonance  fields.  In  the  presence  of  unmistakable  HSZ  effects,  we  expect  to 
see  one  of  the  sets  of  parameters  shown  in  Table  2  emerging  with  unambiguously  non-zero 
values  from  fits  of  the  resonance  field  data.  In  Table  3  we  summarize  the  impact  of  various 
Bn,mUn,m  term  combinations  on  the  parameter  determination;  these  parameter  optimizations 
employed  SH’s  encompassed  by  equation  (5)  and  term  combinations  suggested  by  Table  2.  It  is 
clear  from  these  results  that  HSZ  terms  do  not  significantly  improve  the  quality  of  the  SH 
parameter  fit;  a  small  reduction  in  RMSD  between  observed  and  calculated  resonance  fields  is  to 
be  expected  from  the  addition  of  any  extra  term  to  the  SH.  The  table  shows  that  the  only  term 
with  even  marginal  significance  to  the  fits  is  the  B4, 3X14,3  term  (its  error  being  the  only  one 
smaller  than  the  respective  fitted  parameter  value). 
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We  show  in  Figure  2  the  differences,  calculated  for  both  W-  and  X-band  at  intervals  over 
a  hemisphere  (0  <  0  <  n,  0  <  §  <  u),  between  the  resonance  fields  of  the  3«-»4  AVP  computed 
with  the  (^-independent)  CSH  and  the  corresponding  resonance  fields  computed  with  a  (in¬ 
dependent)  SH  containing  the  HSZ  term  B4>3U4(3.  The  plots  in  Figure  2  (a)  exemplify  the 
differences  computed  for  the  other  five  Cr3*  EPR  transitions  while  the  X-band  differences 
(Figure  2  (b))are  smaller  by  a  factor  of  ~10  and  exhibit  varying  overall  shapes.  The  difference 
map  computed  with  the  B4,-3U4,-3  term  is  related  to  that  shown  for  B4,3U4,3  by  a  shift  of  n/6  along 
the  <f>  coordinate  at  both  fields.  What  may  be  seen  in  Figure  2  (a)  and  Figure  3  is  that  a 
contribution  from  a  term  B4,3U4>3  (with  B4>3  =  0.0005)  creates  deviations  in  a  transition’s  AVP  of 
at  most  +20  G  relative  to  the  CSH  AVP;  this  deviation  is  of  the  order  of  the  RMSD’s  for  SH 
parameter  fits  we  have  calculated  so  far.  Furthermore,  in  consideration  of  Figure  3  (b) 
(especially,  the  flatness  near  0  =  0,n)  and  Figure  3  (a)  (especially,  the  rate  of  change  of  Hres  near 
0=0),  it  is  reasonably  clear  (at  least  for  magnitudes  of  HSZ  parameter  values  encountered  so  far 
in  this  analysis)  that  the  protocol  established  above  for  placing  B0  parallel  to  the  crystal  c  axis  is 
appropriate  for  EPR  spectroscopy  of  the  sapphire:Cr3+  system  at  this  frequency.  In  fact,  several 
calculations  of  GSH  AVP’s  with  several  values  of  B4,.3  s  0.05  show  that  B0  ||  c  always 
corresponds  to  at  least  a  local  extremum  in  the  observed  resonance  field  for  the  3<-*4  transition. 
It  was  precisely  this  concern,  about  the  BS3  effects  on  the  position  of  Cr3+  resonance  field,  which 
led  us  to  initially  describe  the  angles  0  and  (j)  as  “nominal;”  the  angle  assignments  were  based  on 
an  assignment  of  0=0  which  derived  from  the  conventional  z-axis  alignment  approach  described 
above. 

Regarding  possible  errors  in  the  angular  coordinates  at  which  resonance  fields  were 
measured,  a  possible  correction  to  the  data  in  Figure  1  involves  a  systematic  correction  to  the 
value  of  0,  where  it  is  assumed  that  that  the  magnetic  field  is  positioned  at  <j>=0  (the  xz  plane)  for 
all  measurements.  The  functional  forms  for  an  S=3/2  system’s  resonance  field  AVP’s  due  to  the 
sum  of  Zeeman  and  fine  structure  interactions,  in  the  limit  of  high  field  (and  in  the  limit  of  very 

small  higher  order  effects)  are  known  for  <t>=0[31L  A  0  correction  may  be  achieved  by  a  separate 
least-squares  fit  of  the  data  for  each  transition’s  AVP  with  an  adjustable  parameter,  80, 
corresponding  to  a  required  correction  to  the  “nominal”  0  angles.  Analysis  of  the  five  “high 
field”  AVP’s  of  Figure  1  in  this  way  yielded  80  corrections  of  0.00093n,  -0.0005671,  0.0017ti, 
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0.0009571,  and  0.0025tx  (0.17,  -0.10,  0.31,  0.17,  0.45  degrees)  respectively,  and  a  new  angle, 
0=6+80,  is  defined  for  each  resonance  field.  SH  parameter  determination  of  the  0-adjusted 
AVP’s  based  on  equation  (1)  yielded  CSH  parameters  similar  to  those  of  the  previous  analysis, 
but  the  quality  of  the  fit  improved  substantially  (RMSD  =  17.54  G,  see  Table  3). 

It  is  of  interest  to  consider  the  deviations  of  resonance  fields  for  two  calculated  3«->4 
AVP’s  v.hich  differ  only  in  the  value  of  0  assigned  to  each  resonance  field  (80  =  0.3).  The 
differences  between  these  two  sets  of  computed  resonance  fields  as  a  function  of  0  models  the 
deviations  one  would  observe  between  a  resonance  field  data  set  exhibiting  systematic  error  in 
assigning  0  to  resonance  field  measurements  and  a  calculation.  This  is  shown  in  Figure  4.  We 
note  that  this  pattern  of  deviation  bears  some  qualitative  and  quantitative  similarity  to  the 
deviation  in  calculated  resonance  fields  produced  by  including  the  HSZ  term 

(B4  3=0.0005)  in  the  SH  (relative  to  a  calculation  employing  only  the  CSH).  The  latter  deviation 
is  shown  in  Figure  3  (b).  Comparison  of  the  qualitative  features  of  these  deviations  indicates  that 
there  may  be  some  “crosstalk”  between  an  error  in  assigning  0  and  the  effect  of  the  B4,3  HSZ 
interaction  on  resonance  fields  measured  in  the  xz  plane.  Examination  of  Table  3,  however, 
demonstrates  that  0  corrections  did  not  appreciably  change  the  fitted  size  of  HSZ  effects  in  the 
resonance  fields.  The  HSZ  terms  still  contribute  negligibly  to  the  quality  of  SH  parameter 
determination  as  evidenced  by  the  minimal  reduction  in  the  RMSD  between  observed  and 
calculated  resonance  fields  as  HSZ  terms  were  added  to  the  analysis  of  0-corrected  data. 

Farther  angle  corrections  may  be  of  interest  because  we  generally  observed  slight 
resonance  field  differences  in  pairs  of  resonance  field  measurements  that  differ  only  by  a 
nominal  n  rotation  of  the  crystal.  One  of  our  3<-*4  AVP’s  is  exemplary  of  this  phenomenon: 
Bres(0=O,<t>=O)  =  3.0084  T,  Bres(7t,0)  =  3.0092  T.  We  conclude  from  this  that  rotation  of  the  B0 
field  did  not  occur  strictly  in  the  xz  plane,  i.e.,  that  the  crystal  b  axis  (y  direction)  was  not 
parallel  to  the  effective  sample  rotation  axis.  Therefore  some  (additional)  correction  to  the 
nominal  angle  assignments  may  be  appropriate.  Small  angle  corrections  addressing  this  problem 
may  be  computed  by  noting  that  the  resonance  fields  of  the  outer  AVP’s  are  proportional  to 
cos2(0)  for  0  values  near  0  and  n  (in  the  limit  of  very  small  higher  order  interactions): 

Bres(0)=^+5COS2(0)  (7) 
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The  0=O,7i  resonance  fields  may  be  used,  in  conjunction  with  a  cos2(0)  fit  of  the  whole 
3<-»4  a  VP,  to  compute  an  apparent  0  for  a  ti  rotation  about  the  effective  rotation  axis  (the  actual 
value  of  0  and  the  rotation  angle,  7t,  are  different  because  the  rotation  axis  is  not  exactly  parallel 
to  the  crystal  b  axis).  By  assuming  that  resonance  fields  near  the  yz  plane  vary  as  cos  (0),  we 
have: 

Bres(0  rotation)  —  Bres(7T  rotation)  =  0.0008  Tesla  =  i?(cos2(0)  -  cos2(0n))  (8) 

where  0,  is  the  angle  B0  makes  with  the  z-axis  after  a  n  rotation  about  the  effective  rotation  axis, 
and 

cos(0„)  =  -|  [(0.0008/2?)  +  1],/2 1  (9) 

where  the  value  of  B  was  obtained  to  a  good  approximation  from  the  cos2(0)  fit  of  the 
whole  3<-*4  AVP.  Using  the  actual  rotation  angle  (tt)  and  the  effective  0n  computed  in  this  way 
we  may  employ  a  rotation  matrix,  which  depends  only  on  the  effective  rotation  axis  vector 
components  (referred  to  the  laboratory  frame)  and  the  actual  angle  of  rotation  about  the  effective 
axis,  to  compute  two  components  of  the  rotation  axis  vector.  We  make  the  simplifying 
assumption  that  the  rotation  axis  lies  in  the  yz  plane.  With  knowledge  of  the  effective  rotation 
axis  we  may  discern  the  corrected  values  of  0  and  4>  for  resonance  fields  of  each  rotation  position 
in  the  AVP.  This  angle  correction  -  applied  only  to  the  one  3+->4  AVP  mentioned  above  as  a 
supplement  to  the  0  +  80  corrections  mentioned  above  (applied  previously  to  all  AVP’s)  - 
improved  the  overall  SH  parameter  fit,  giving  a  RMSD  of  16.00  G  (see  table  Table  3). 

We  note  that  since  this  angle  correction  positions  some  data  away  from  the  xz  plane  (to 
which  data  were  nominally  assigned  at  first),  the  determination  of  the  B4,-3  parameter  value  is 
now  not  accompanied  by  a  large  error.  Its  value  (from  the  EPRNMR  parameter  optimization) 
increased  dramatically  as  well.  The  values  of  the  CSH  parameters,  however,  were  not 
appreciably  changed,  indicating  their  robustness  vis  a  vis  our  data,  the  several  HSZ  scenarios 
considered,  and  this  0  and  <j)  angle  correction  employed.  This  robustness  may  be  good  evidence 
for  the  (^-independence  of  the  EPR  spectra.  If  the  spectra  are  ^-independent,  detectable  HSZ 
effects  associated  with  SH  terms  B4+3U4+3  should  be  absent.  Note  that  effects  due  to  terms 
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B2,oU2,o  and  B4>0U4,o  are  inherently  (^-independent  and  may  therefore  not  be  ruled  out  by  such 
geometric  considerations. 

In  addition  to  producing  reasonable  values  of  the  CSH  parameters,  the  analysis  to  this 
point  serves  to  underscore  the  importance  of  including  data  that  falls  outside  principal  planes  in 
the  assessment  of  HSZ  parameters.  When  data  is  not  obtained  outside  of  the  principal  planes, 
some  dimensions  of  the  SH  parameter  space  are  left  unexplored  due  to  the  periodic  zero 
crossings  -  some  of  which  lie  in  principal  planes  -  of  trie  HSZ  term  contributions.  Examples  of 
such  zero-crossings  may  be  seen  in  Figure  2  (a).  With  the  goal  of  gaining  a  more  secure 
foundation  on  which  to  include  or  neglect  HSZ  effects  in  the  interpretation  of  the  EPR  of  the 
sapphire:Cr3+  system,  we  now  consider  some  resonance  field  data  collected  outside  of  the  xz 
plane. 

With  the  present  sample  and  sample-rotation  apparatus,  we  are  unable  to  effect  full 
rotations  of  B0  in  planes  far  from  the  plane  perpendicular  to  the  fiber  axis.  Nevertheless,  we  can 
rotate  B0  away  from  this  plane,  by  as  much  as  rbt/9,  by  means  of  the  microwave  coupling 
adjustment.  These  rotations  are  fruitful  because  the  signals  of  interest  are  of  more  than  adequate 
strength  even  when  microwave  coupling  is  suboptimal,  and  because  a  rotation  of  this  extent  does 
not  substantially  decrease  the  component  of  Bi  ±  B0.  Resonance  fields  were  measured  at  various 
coupling  adjustments  (away  from  the  position  where  B0  is  nominally  perpendicular  to  the 
sample’s  long  axis,  i.e.,  at  positions  where  B0  is  nominally  out  of  the  xz  plane)  for  many 
rotational  positions  of  the  sample  about  its  long  axis.  For  such  measurements,  the  (0,<)>) 
orientation  of  B0  in  the  crystal  may  be  calculated  from  the  cumulative  effects  of  rotations  of  the 
crystal  about  two  axes  in  the  laboratory  reference  frame:  the  sample’s  cylinder  axis  (oriented 
initially  along  the  laboratory  frame  y-direction)  and  the  coupling  rotation  axis  (oriented  along  the 
laboratory  frame  x-axis).  We  measured  resonance  fields  (as  many  as  three  of  the  “strong”  Cr3+ 
EPR  transitions)  at  over  100  field  positions  in  the  crystal  frame,  and  these  positions  are  shown  in 
Figure  5. 

The  results  of  spin-Hamiltonian  parameter  optimization  on  this  set  of  297  in-  and  out-of- 
plane  resonance  fields  are  listed  in  Table  3.  An  initial  optimization  was  performed  employing 
the  CSH  and  no  angle  corrections,  and  the  overall  quality  of  the  parameter  optimization  is 
indicated  in  the  table  (RMSD  »  39.85  for  297  data  points).  This  optimization  clearly  represents 
an  inferior  parameter  optimization  vis  a  vis  even  the  optimizations  for  angle-uncorrected  in-plane 
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date  reported  in  Table  3.  This  may  indicate  the  need  for  resonance  field  orientation  corrections. 
Further  support  for  this  hypothesis  is  indicated  by  the  dramatic  improvement  achieved  in  a 
parameter  fit  minimization  in  which  all  elements  of  the  g  and  D  matrices  were  allowed  to  vary 
(RMSD  =  22.39).  Several  orientation  corrections  to  this  data  are  now  considered. 

Parametric  variation  of  0a  (the  raw  azimuthal  angle,  i.e.,  the  angle  of  rotation  about  the 
sample’s  cylinder  from  the  point  where  B0  ||  z)  and  y  (the  raw  coupling  angle),  the  angles  used  to 
calculate  0  and  <|>  direction  of  B„  in  the  crystal,  reveal  only  small  improvements  to  fits  when  new 
angles  ©  and  O  are  defined  (derived  from  0a  — ►  0a  +  80a  +  £e  *0a  and  y  — y  +  8y  +  £y*y,  where 
parameters  were  varied  in  these  ranges:  -1.0  <  80a,  8y  <  +1.0  and,  -0.06  <  £e,  £y  <  0.06).  The 
angle  corrections  which  produced  both  the  minimum  RMSD  and  the  minimum  rhombicity  of  the 
D  tensor  (under  optimization  in  which  all  D  matrix  elements  are  allowed  to  vary)  are:  80a  = 
0.75, 8y  =  0.0,  se  =  0.0,  &y  =  0.03. 

Our  CSH  parameter  values  are  in  good  agreement  with  the  values  reported  by  Chang,  et 
al.  [12]  We  have  noted  that  these  values  are  stable  when  fits  of  HSZ  terms  are  attempted  and 
suggested  that  this  may  be  due  to  the  lack  of  significant  HSZ  effects  in  the  resonance  fields.  The 
quantum  mechanical  perturbation  treatment  of  electronic  states  shows  that  HSZ  effects  trilinear 
in  electronic  spin  vector  components  are  proportional  to  (X/A)3  (where  X  is  the  spin  orbit 
coupling  constant  for  Cr3+  and  A  is  the  crystal  field  splitting).  In  consideration  of  published 

literature  in  which  HSZ  effects  are  observed!^],  and  in  consideration  of  the  quantum 
mechanical  origin  of  these  terms,  it  is  likely  that  paramagnetic  ions  with  ever,  marginally  higher 
spin-orbit  couplings  (X)  are  better  candidates  than  Cr3+  for  observation  of  HSZ  effects. 

We  have  alluded  to  previous  measurements  of  HSZ  effect st^  2]  and  wish  to  further 
comment  on  the  relative  success  these  authors  enjoyed.  The  sapphire:©!-3*  system  possesses 
electronic  quadrupole  splitting  which  is  evidenced  by  the  substantial  variation  of  the  resonance 
fields  of  the  outer  (3/2<-»l/2,  -l/2<->-3/2)  transitions  as  a  function  of  crystal  orientation  in  the 
magnetic  field.  But  we  have  also  noted  some  similarity  between  the  errors  in  the  AVP’s  that  are 
produced  by  angular  misalignment  and  the  computed  effects  in  the  AVP’s  due  to  addition  of 
certain  HSZ  terms.  The  work  on  Co  shows  that,  with  D=0  (rigorously),  the  HSZ  effect  appears 
on  an  otherwise  stationary  Bres  background  and  is  thus  easily  discernible  without  consideration  of 
small  angular  location  errors.  On  the  other  hand,  the  determination  of  small  HSZ  effects  from 
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the  Cr3+  system  appears  to  require  forbidding  angular  precision  given  the  large  variation,  due  to 
the  fine  structure  interaction,  of  resonance  fields  with  angle. 

The  relationship  between  the  values  of  the  readily  observed  EPR  parameters  (g,,  gj.,  D) 
and  the  arrangement  of  the  electronic  states  of  Cr3+  in  a-Al203  has  been  discussed  thoroughly  by 
Macfarlane[33-35]  jn  the  context  of  a  crystal  field  model.  It  was  already  apparent  from 
published  experimental  data  that  the  “textbook”  relationship  (D  =  A.(g,  -  gj.),  X  =  spin-orbit 
coupling  parameter)  was  not  applicable  to  Cr3+  in  sapphire.  It  was  with  Macfarlane’s  work  and 

that  of  Sugano  and  TanabePG]  that  the  complexity  of  the  quantum  mechanical  relationship 
between  D  and  the  g-values  began  to  be  rationalized.  This  complexity  arises  from  the  fact  that 
the  second-order  contributions  that  lead  to  the  textbook  expression  vanish  in  sapphire:Cr3+,  so 
that  only  in  higher  order  calculations  are  the  origin  of  and  relationship  between  D  and  g-values 
revealed.  Exact  calculations  of  D  and  g  are  extremely  difficult;  the  theoretical  derivations 
provided  in  the  work  of  Macfarlane  show  that  D  and  g  arise  from  numerous  small  high-order 
contributions.  Even  so,  the  optical  spectrum  provides  X  and  the  Racah  parameters  B,  and  C  so 
that  decent  approximations  of  the  EPR  parameters  may  be  found.  We  note  in  this  context  that 
Rimmer  and  Johnston[37]  identified  small  adjustments  to  the  parameter  values  of  Macfarlane’s 
model  that  allowed  much  improved  agreement  of  the  computations  of  EPR  parameters  with 
experimentally  observed  values. 

In  view  of  the  complexity  of  the  likely  models  of  Cr3+  in  a  trigonal  field,  Macfarlane, 
commenting  on  the  analysis  of  Tanabe  and  Sugano,  pointed  out  the  disagreeable  nature  of 
including,  ad  hoc ,  spin-orbit  coupling  anisotropy  and  multiple  spin-orbit  coupling  parameters  in 
the  treatment  of  this  system.  Even  so,  Klein,  et  al. ,  introduced  these  phenomenological 
parameters  and  was  able  to  rationalize  the  (probably  incorrect)  g-values  of  Manenkov  and 

Prokhorov!^*  2],  We  find,  however,  that  if  we  attempt  to  analyze  the  degree  of  spin-orbit 
coupling  anisotropy  using  Klein’s  model  (which  posits  spin-orbit  coupling  anisotropy)  and  our 
(better)  g-values  which  exhibit  small  Ag,  we  calculate  the  spin-orbit  anisotropy  to  be  0.5%  (a 
factor  of  32  smaller  than  the  figure  of  Klein).  In  other  words,  our  measurements  of  Ag,  in 
agreement  of  those  of  several  other  groups,  predict  essentially  isotropic  spin-orbit  coupling 
within  the  Klein  model.  Macfarlane’s  simpler  analysis  is  therefore  to  be  preferred  as  it  finds  the 
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correct  g-anisotropy.  This  position  is  further  sustained  if  we  consider  this  model  in  light  of  the 
better  EPR  parameter  values,  e.g.,  those  of  Chang,  et  all1 2]  or  those  of  the  present  work 

Conclusion 

We  have  presented  the  results  of  a  study  of  the  orientation  dependence  of  94.9  GHz  Cr3+ 
EPR  cr.  ectra  in  a  single-crystal  of  neutron-irradiated  white  sapphire.  We  discussed  the  results  of 
previous  EPR  studies  of  this  system,  and  the  heretofore  universal  application  of  the  CSH  to  the 
interpretation  of  these  results.  We  examined  a  key  weakness  of  the  CSH,  the  absence  of  HSZ 
terms  which  are  in  principle  relevant  on  group  theoretical  grounds  to  the  study  of  S  =  3/2 
systems.  The  prospects  for  the  high-field  EPR  determination  of  the  size  of  HSZ  effects  in  the 
sapphire:Cr3+  system  were  analyzed,  and  we  discussed  the  significant  benefit  that  positive 
identification  of  such  effects  would  entail,  namely,  the  more  complete  elucidation  of  the  point 
symmetry  of  the  paramagnetic  ion  in  the  crystal.  This  investigation  revealed  no  conclusive 
evidence  of  HSZ  effects;  in  light  of  this  result,  we  discussed  possible  alterations  and 
improvements  of  experimental  conditions  that  may  facilitate  the  identification  of  HSZ  effects, 
namely,  improvement  of  angular  measurements  and  the  use  of  measurements  out  of  the  xz  plane. 
We  believe  that  unambiguous  measurements  of  HSZ  effects,  if  fruitful,  would  constitute  a 
significant  achievement  for  high-field  EPR  spectroscopy.  The  ability  to  measure  HSZ  effects 
could  have  much  to  offer  in  the  analysis  of  chemical  structure  problems  from  biology  to  material 
sclent. 
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Figure  Captions 


Figure  1  Angular  variation  of  the  three  “high  field”  94.9  GHz  EPR  resonance  fields  of  Cr3 
in  nominally  pure  neutron-irradiated  thermally  annealed  01-AI2O3  (sapphire)  single  crystal. 
Rotation  occurred  (nominally)  about  the  crystal  b  axis  (<j>=0);  at  0=0,  Bext  is  parallel  to  the  crystal 

c  axis.  The  five  tracks  include  three  for  the  3+-*4  transition  (0,13,0,  lowest  Bres  at  0=0),  and  one 
each  for  the  2<-*3  transition  (A,  intermediate  Bres  at  0=0)  and  the  l«-»2  transition  (v,  lowest  Bres 
at  0=7i/2). 

Figure  2  Calculated  angular  variation  (a)  at  W-band  (94.9  GHz)  and  (b)  at  X-band  (9.49 
GHz)  of  the  contribution  of  the  B4)3U4,3  SH  term  to  the  resonance  field  of  the  3«-»4  transition  of 
an  S=3/2  system  with  SH  parameters  gj_  =  1.9814,  gj  =  1.9810,  D  =  -5.7413  GHz,  B4,3  =  0.0005 

(dimensionless!^]). 

Figure  3  The  calculated  angular  variation  pattern  of  the  3«-*4  transition  resonance  fields  (0 
s  ©  <  7i, <(>=0)  of  (a)  the  CSH;  (b)  the  slice  of  the  difference  plot  of  Figure  2  (a)  for  <|>  =  0;  (c)  the 
slice  of  the  difference  plot  of  Figure  2  (a)  (0  =  2rc/3,.  0  £  <j>  s  tc).  Figures  (b)  and  (c)  show  the 
maximum  possible  contributions  of  the  B4,3U4>3  term  to  the  ((^-independent)  Bres  surface 
exemplified  by  the  (<|>=0)  slice  (a). 

Figure  4  The  difference  between  two  3<-*4  CSH  AVP’s  computed  with  gx  =  1.981 1,  gj  = 
1.9820,  D  =  -5.7413  GHz.  AVP’s  differ  in  that  one  is  shifted  along  the  0  coordinate 
( 0-»0+O.OO1677t).  The  plot  indicates  the  effect  of  a  systematic  error  in  assigning  0  to  sample 
positions  and  should  be  compared  to  Figure  3  (b). 

Figure  5  The  orientations,  in  the  crystal’s  reference  frame,  of  the  applied  magnetic  field  for 
294  rotation/coupling  variations. 
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Tables 


V  (GHz) 

Sii 

gx 

D  (GHz) 

Reference 

8.96,  11.97 

1.984±0.0006 

1.9867±0.0006 

5.732*f 

[1,2] 

37.86 

1.982±0.002 

1.979±0.009 

-5.732±0.03 

[3,4] 

9.309,  11.593 

2.003±0.006 

2.00±0.02 

5.79±0.03 

[5] 

9.58 

1.98* 

1.98* 

5.7*t 

[6,  7] 

9-18 

~1.982*J 

~1.982*J 

5.74§ 

[9] 

24.2 

1.97*J 

1.97*| 

5.807*f 

[38] 

9.5,  24 

1.9817±0.0004 

1.9819±0.b006 

-5.7465±0.002 

[12] 

94.9 

1.9812xxx 

1.9820xxx 

5.737xxx 

This  work 

Table  1  Room  temperature  EPR  studies  sapphire:Cr3+,  with  reported  spectrometer  frequency  and 
spin  hamiltonian  parameter  measurements. 

f  value  reported  as  absolute  value 

*  no  error  reported 

$  only  one  g-value  reported 

§  value  extracted  from  plot  figure 


Symmetry 

Compliant 

Spin 

Hamiltonain 

terms 

TC=H: 

+B2,0 

+B4;0 

+B2,0 

+B4,o 

+B4;±3 

+B2,0 

+B4)o 

+B4)+3 

+B4j.3 

Point  group 

C  3h 

Civ 

c3 

symmetry 

D 3h , 

D  3 

s6 

D  3d 

Table  2  The  SH  terms  from  equation  (10-4)  which  comply  with  the  symmetry  restrictions  of 
several  point  groups. 
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matrices _ 

Fuliy  1.9808  19817  1.9821  5.705  ;£  0.0077  9.4e-6  -2.5e-5  -6.9e-5  -8.1e-5  0  0  23.65  294 

unconstrained  ±0.0006  ±0.0002  ±0.0008  ±0.009  ±6.9e-5  ±7.2e-5  ±7.7e-5  ±7.7e-5 


angle  offset  text)  text) 

optimization _ _ _ _ _ _ _ 


Table  3  Root  mean-square  deviations  (RMSD)  of  simulated  resonance  fields  from  observed 
resonance  fields.  Bi,m  are  HSZ  parameters  discussed  in  the  text.  See  the  text  for  information 
about  50  values  used  in  certain  fits.  Extremely  large  errors  ( e.g .,  data  row  3,  data  col.  9)  indicate 
that  the  impact  of  that  parameter  on  the  computed  resonance  fields  is  negligible  for  the  angular 
orientations  of  the  fitted  data. 


III-20 


References 


[1]  A.  A.  Manenkov,  A.  M.  Prokhorov,  J.  Exp.  Theor.  Phys.  (USSR)  28  (1955)  762. 

[2]  A.  A.  Manenkov,  A.  M.  Prokhorov,  Sov.  Phys.  JTEP  1  (1955)  611. 

[3]  G.  M.  Zverev,  A.  M.  Prokhorov,  J.  Exp .  Theor .  Phys.  (USSR)  34  (1958)  513. 

[4]  G.  M.  Zverev,  A.  M.  Prokhorov,  Sov .  Phys.  JTEP  7  (1958)  354. 

[5]  J.  E.  Geusic,  Phys.  Rev.  102  (1956)  1252. 

[6]  M.  M.  Zaripov,  1. 1.  Shamonin,  J.  Exp.  Theor.  Phys.  (USSR)  30  (1956)  291. 

[7]  M.  M.  Zaripov,  1. 1.  Shamonin,  Sov.  Phys.  JTEP  3  (1956)  171. 

[8]  E.  O.  S.-D.  Bois,  Bell  Sys.  Tech.  J.  38  (1959)  271. 

[9]  H.  Klein,  U.  Scherz,  M.  Schulz,  H.  Setyono,  K.  Wisnewski,  Z.  Phys.  B  28  (1977)  149. 

[10]  G.  Kido,  N.  Miura,  Appl  Phys.  Lett.  41  (1982)  569. 

[1 1]  R.  W.  Terhune,  J.  Lambe,  C.  Kikuchi,  J.  Baker,  Phys.  Rev.  123  (1961)  1265. 

[12]  T.  Chang,  D.  Foster,  A.  H.  Kahn,  J.  Res.  NBS  83  (1978)  133. 

[13]  A.  Abragam,  M.  H.  L.  Pryce,  Proc.  Royal  Soc.  A205  (1951)  135. 

[14]  A.  B.  Roitsin,  Phys.  Stat.  Sol  (b)  104  (1981)  11. 

[15]  D.  G.  McGavin,  W.  C.  Tennant,  J.  A.  Weil,  J.  Mag.  Res.  87  (1990)  92. 

[16]  H.  A.  Buckmaster,  R.  Chaterjee,  Y.  H.  Shing,  J.  Mag.  Res.  4  (1971)  85. 

[17]  V.  G.  Grachev,  Sov.  Phys.  JTEP  65  (1987)  1029. 

[18]  S.  Geschwind,  Phys.  Rev.  121  (1961)  363. 

[19]  G.  F.  Koster,  H.  Statz,  Phys.  Rev.  113  (1959)  445. 

[20]  H.  A.  Buckmaster,  R.  Chaterjee,  Y.  H.  Shing,  Phys.  Stat.  Sol.  A  13  (1972)  9. 

[21]  B.  Bleaney,  Hyperfine  Interactions ,  Academic  Press,  New  York  1967. 

[22]  B.  Bleany,  Proc.  Phys.  Soc.  A73  (1959)  939. 

[23]  S.  A.  Alt'shuler,  B.  M.  Kozirev,  Electron  Paramagnetic  Resonance  in  Compounds  of  Transition  Elements , 
John  Wiley  &  Sons,  New  York  1972. 

[24]  W.  J.  C.  Grant,  M.  W.  P.  Strandberg,  J.  Phys.  Chem.  Sol.  25  (1964)  635. 

[25]  G.  Bacquet,  J.  Dugas,  C.  Escribe,  J.  M.  Gaite,  J.  Michoulier,  J.  Chem.  Phys.  C7  (1974)  1551. 

[26]  J.  Michoulier,  J.  M.  Gaite,  J.  Chem.  Phys.  56  (1972)  5205. 

[27]  D.  S.  Schonland,  Molecular  Symmetry:  An  Introduction  to  Group  Theory  and  its  Uses  in  Chemistry ,  Van 
Nostrand  Reinhold  Co.,  London  1965. 

[28]  W.  M.  Yen,  Preparation  of  Single  Crystal  Fibers ,  World  Scientific,  Singapore  1995. 

[29]  D.  A.  Schwartz,  E.  D.  Walter,  S.  J.  Mcllwain,  V.  N.  Krymov,  D.  J.  Singel,  Appl.  Magn.  Reson.  16  (1999) 
223. 

[30]  M.  J.  Mombourquette,  J.  A  Weil,  D.  G.  McGavin, ,  University  of  Saskatchewan,  Saskatoon  1993. 

[31]  A.  Abragam,  B.  Blear. zy,  Electron  Paramagnetic  Resonance  of  Transition  Ions ,  Clarendon  Press,  Oxford 
1970. 

[32]  F.  S.  Ham,  G.  W.  Ludwig,  G.  D.  Watkins,  H.  H.  Woodbury,  Phys.  Rev.  Lett.  5  (1960)  468. 

[33]  R.  M.  Macfarlane,  J.  Chem.  Phys.  39  (1963)  3118. 

[34]  R.  M.  Macfarlane,  J.  Chem.  Phys.  42  (1965)  442. 

[35]  R.  M.  Macfarlane,  Phys.  Rev.  B  1  (1970)  989. 

[36]  S.  Sugano,  Y.  Tanabe,  J.  Phys.  Soc.  Jap.  13  (1958)  880. 

[37]  D.  E.  Rimmer,  D.  F.  Johnston,  Proc.  Phys.  Soc.  89  (1966)  943. 

[38]  C.  Kikuchi,  J.  Lambe,  G.  Makhov,  R.  W.  Terhune,  J.  Appl.  Phys.  30  (1959)  1061. 


III-21 


Appendix  IV 


C++  Code  for  Simulation  of  Rapid  Passage  Effects  in  EPR 


Numerical  Simulation  of  the  Bloch  Equations 

In 

CW-EPR 

I.  Introduction 

This  document  discusses  methods  to  simulate  CW-EPR  spectra  through  numerical 
integration  of  the  bloch  equations. 

II.  Equations  &  Theory 

The  magnetization  vector,  M ,  is  defined  as: 

(K) 

M„ 


When  a  collection  of  spins  resting  at  equilibrium  M°  =  0  are  submitted  to  a 

Ik) 

modulated  magnetic  field  pit)  (shown  below)  along  the  z-axis  and  microwaves  at  amplitude 
power  Hi  along  the  x-axis  of  the  rotating  frame,  their  behavior  is  defined  via  the  following 
differential  equations: 

(  1  ^ 

~  ~r-p(0  0  f  1 

T2  0 

Y'Pif)  -Y'h\  -M+  o 


y-Hx 


pit)  =  H0-  (Hc  +  Hm  cos(zzrm  •  0) 

Where  T,  and  T2  are  the  relaxation  rates,  y  is  the  gyro-magnetic  ratio  for  a  free  electron. 
Ho  is  the  resonant  field,  He  is  the  current  field,  com  is  the  modulation  frequency  (rad/sec),  and  Hm 
is  the  modulation  amplitude.  The  Jacobian  is  obtained  below  for  use  with  other  simulation 
methods,  such  as  LSODES,  a  stiff  numerical  integration  method  in  the  ODEPACK  bunch  of 
FORTRAN  programs. 


dMx 
dt  ■  dMx 
dMy 
dt  ■  dMx 
dMz 
8t-8Mr 


dMx 

8MX 

dt  ■  8My 

dt  ■  8MZ 

8My 

8My 

dt  ■  8My 

dt  •  8MZ 

8MZ 

8MZ 

dt  •  8My 

dt  ■  dM. 
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J= 


r,m  8M> 

T2  3Mz  r  3MZ  y  3MZ 


I  dMz  ^  dp(t)  M 


T2  3My  3My 


y 

drp(t).,  , ,  BM.  \  SMt  dM,  8yp(t)  , ,  3MZ  1  3My 

^■^nrrTjd-r^'iir,  isr"-+w<08«,  r,«. 


dMy 


3My 


1  3Mz 


3MZ  T,  3MZ 


3My  1  dMz 
3My  ~  fl  3My 


y-rrt)^ 

T2  3M.  3Mz  y  ■ 


3My 

'■m. 


|£- 

dMy  3MZ  dM. 


y»x 


dMv 


I  3My 
'  T2  3Mz 
1  dMz 


■ m . 


dM.  7]  3MZ 


dM. 

dM. 


f 


J  = 


__L 

t2 

yp(t) 

o 


-yp(t)  0 


1 


- yH \ 
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The  line-shape  function  L(H)  will  be  defined  as  (see  above): 


-In  2 


H„-H 


L{H)  =  e 

Others  line-shape  functions  may  later  be  implemented  (e.g.  Lorentzian,  Voigt)  and  or 
powder  patterns  generated  by  program  and/or  experiments.  This  document  will  deal  with  the 
basic  gaussian  function  so  that  the  implementation  can  be  checked  for  errors  and/or 
inconsistencies  before  writing  additional  functionality. 


The  solutions  to  the  equations  at  different  offset  values  (He)  are  called  wavepackets  or 
spinpackets  throughout  this  document.  A  multitude  of  packets  with  He  ranging  from  -H^  to 
+Hmax  in  increments  of  HI  (G)  or  \/T2  (G),  depending  on  which  is  smaller  (increment  va.ue  is 
deemed  dH),  are  created  to  be  used  to  simulate  a  spectrum  using  the  equations  e  ow. 
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Hi  is  then  defined  as 

Hi  =  Hs  j  -Hcj 

The  spinpackets  are  then  summed  up  over  the  line  shape  function  through  the  equation: 

n-number  of  _  packets 

( 1 )  suntj  (t)  =  !(//,)•  packets ,  (t)-dH 

i=0 

The  spectrum  for  the  0  &  90  degrees  out  of  phase  is  then  extracted  by  calculating  the 
fourier  coefficients  for  each  pt(j)  on  the  line  using  (2)  &  (3) 

^ma\  2  *  dt 

(2)  spectra**  j  =  YjSumft)  •  cos(/z  -&J)-  j — — r 

t=t0  Vmax  ‘o/ 

^max  2  *  dt 

(3)  spectra90* j  =  •  cos(/i  •  tzrmr  +  -)  •  j — — -r 

^  Vmax  / 

■  *« 

h>\ 

Spectra0,h  and  spectra90'11  is  calculated  for  the  M*  and  My  packet  set.  Small  h  m  these 
equations  is  the  integral  harmonic  spectrum  (relative  to  modulation),  to  be  extracted. 

III.  Programming  Method 

Using  fourth  order  Runge-Kutta  integration  methods,  a  program  was  written  in  Matlab  to 
produce  a  plot  of  one  Wavepacket,  Mx,My,  and  Mz  as  a  function  of  time.  Later,  an  offsets  were 
added  for  the  acquisition  of  an  evenly  spaced  array  of  wavepackets,  with  the  offset  depending  on 
HI  (G)  or  1/T2  (G),  whichever  was  smallest.  Due  to  the  relative  inefficient  of  the  interpreted 
language  in  Matlab,  a  C  version  of  the  code  on  a  UNIX  system  was  implemented,  thus  giving  rise 
to  the  blochc  programs.  As  a  result,  wavepackets  were  simulated  in  ~360  fold  less  time  than  the 
Matlab  product.  Further  improvements  were  added  to  automatically  catch  and  handle  when  the 
integration  “blows-up”,  i.e.  the  magnetization  vector  components  absolute  values  become  larger 
than  a  programmed  ceiling  value. 

After  the  method  for  simulating  multiple  wavepackets  in  a  reasonable  amount  of  time 
-;,  as  formulated,  a  new  program  called  spec  was  implemented.  This  program  takes  the  output 
wavepackets  from  the  blochc  program.  After  subtracting  the  packet  offsets,  for  each  point  on  the 
spectrum,  spec  sums  the  wavepackets  with  the  gaussian  function  weight  centered  around  the 
selected  point  and  demodulates  the  composite  Wavepacket  time  to  obtain  the  Fourier  coefficient. 
The  Fourier  coefficients  are  calculated  for  each  point  on  the  spectral  line.  Using  the 
demodulation  function,  a  nth  harmonic  spectrum  can  also  be  produced.  The  current  version  of 
spec  only  implements  the  1st  and  2nd  harmonic  of  the  calculation. 

IV.  Performance  &  Data  Analysis 

To  quantify  the  performance  of  the  blochc  program,  several  factors  have  to  be 
considered.  Some  these  factors  include,  but  are  not  limited  to,  modulation  frequency,  modulation 
amplitude,  T2  values,  and  Wavepacket  offset.  These  all  effect  the  step  size  required  by  the 
Runge-Kutta  method  to  produce  a  “reliable”  result.  This  section  will  play  with  the  parameters  of 
the  blochc  program  while  measuring  the  amount  of  time  required  to  finish  the  Wavepacket.  First, 
we  will  deal  with  a  spinpacket  with  offset  of  zero,  vaiying  the  parameters  that  may  affect  the 
performance  of  the  program.  Nine  runs  will  be  averaged  to  obtain  the  proper  spectrum.  The 
computer  used  will  be  an  Octane  250MHz  computer  running  DEC  UNIX  (version?).  The 
parameters  and  results  for  the  center  wavepackets  are  shown  below  in  Table  I  below. 
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Table  I  -  Performance  Analysis  Using  the  Central  Spinpacket 


Run  # 

HI  (G) 

T1 

(sec) 

T2 

(sec) 

Hv 

(Hz) 

Hm  (G) 

Cycles 

Date/Time 

Time 

Completed 

1 

0.065 

18e-6 

le-6 

500 

0.5 

3 

2 

0.065 

18e-6 

le-6 

1000 

0.5 

3 

3 

0.005 

18e-6 

le-6 

500 

0.5 

3 

4 

18e-6 

le-6 

1000 

0.5 

3 

5 

3  6e-6 

le-6 

500 

0.5 

3 

6 

0.065 

36e-6 

le-6 

1000 

0.5 

3 

7 

0.005 

36e-6 

ler6 

500 

0.5 

3 

8 

0.005 

36e-6 

le-6 

1000 

0.5 

3 

9 

0.065 

18e-6 

le-6 

500 

1.0 

3 

10 

0.065 

18e-6 

le-6 

1000 

1.0 

3 

11 

0.005 

18e-6 

le-6 

500 

1.0 

3 

12 

0.005 

18e-6 

le-6 

1000 

1.0 

3 

IV.  Conclusion 

A.  Usage  of  the  blochc4  program 

This  is  the  command  line  for  running  blochc4. 

>  Blochc4  <packets>  <substeps>  <Jpts>  <Hm>  <v>  <H1>  <T1>  <T2>  <cycles>  <filename> 


<packets> 

(integer)  The  number  of  packets  to  calculate.  The  packets  are  offset  by  HI  at  each 
iteration,  i.e.  Offset  of  1st  packet  will  be  (H1(G)  or  l/T2(G))*(packets/2-l),  center 
packet  will  have  an  resonance  offset  of  zero.  This  will  only  apply  if  <packets>  is 
an  odd  integer,  otherwise,  no  calculation  of  the  zero  offset  will  be  performed. 

<substeps> 

(integer)  Initial  number  of  substeps  between  points  to  take.  Program  will  adjust 
this  number  to  obtain  a  “non-blowup”  result  of  the  simulation. 

<fpts> 

(integer''  Mumber  of  Doints  to  save  to  a  file.  File  format  will  be  discussed  below. 

<Hm> 

(real)  Y>  the  desired  modulation  height  (G)  for  the  simulation.  Okay  to  set  to  zero 
if  no  modulation  is  desired. 

<v> 

(real)  Modulation  frequency  (Hz) 

<H1> 

(real)  Microwave  power  amplitude  (G). 

<T1> 

(real)  Relaxation  Rate  (sec) 

<T2> 

(real)  Relaxation  Rate  (sec) 

<cycles> 

(integer)  Number  of  modulation  cycles  to  run  simulation.  Must  be  at  least  one. 

1 - 

<filename> 

(flic  string)  filename  of  where  the  results  are  to  be  saved.  When  blochc4  is  done, 
the  results  are  saved  in  the  ./work  directory  of  where  blochc4  is  run.  For  proper 
operation,  make  sure  that  the  ./work  directory  exists. 

- l - 1 - - - — — : .  1  TTTt  • n  _ L  1 T yl 

rue  iormai  oi  Diocncn  prugiam  m  ™ - - - 

three  files.  Two  of  these  files  are  <filename>x.red  and  <filename>y.red  which  are  the  Mx  and 
My  components  of  the  magnetization  vector  as  a  function  of  time.  In  the  first  row  of  the  file  is 
the  offsets  that  were  used  for  each  Wavepacket.  The  1st  column  contains  the  time  vector  that  was 
used.  The  remaining  columns  are  for  the  magnetization  components  that  were  calculated  for  each 
Wavepacket.  The  third  file  contains  the  parameters  that  were  used  to  calculate  the  simulation 
(<filename>.par). 


B.  Usage  of  spec  programs 

This  is  the  command  line  for  running  the  spec  program. 

>  spec  <basefile>  <points>  <linewidth  >  <bfield>  <efield>  <Center> 


<basefile> 

(filename  string)  Wavepacket  filename.  Notice  that  this  is  the  basefilename,  i.e., 
when  calculating  a  spectrum  from  the  wavepackets  generated  using  <filename> 
above,  only  use  the  filename  as  entered  into  the  blochc4  program.  This  program 
will  append  to  appropriate  extension  and  load  in  the  Mx  and  My  components  of 
the  simulation.  Also  note  that  spec  expects  to  read  from  the  same  ./work  directory 
as  blochc4  does. 

<points> 

(integer)  Number  of  data  that  you  want  calculated  with  the  spectrum.  (1000 
usually  used).  _ 

<linewidth> 

(real)  Linewidth  of  the  gaussian  function  (1/2  height)  of  which  you  are  calculating 

(real)  beginning  field  of  spectrum.  Note  :  this  needs  to  be  less  than  <efield>. 

(real)  ending  field  of  spectrum. 

<Center> 

(real)  point  in  the  spectrum  where  the  gaussian  function  reaches  its  maximum 
(central)  value. 

This  program  will  generate  two  files,  one  for  the  spectra  calculated  using  the  Mx 
component  (  called  <basefile>x.spec)  and  the  other  for  the  My  component  (called 
<basefile>y.spec). 


C.  Usage  of  bspec  scripting  program 
Here  is  the  command  line  for  the  bspec  UNIX  script. 

Bspec2  <packets>  <subpts>  <packetpts>  <Hm>  <v>  <H1>  <T1>  <T2>  <cycles> 
<spec pts>  <linewidth>  <bfield>  <efxe\d>  <center>  <basefile> 

Bspec2  combines  the  functionality  of  blochc4  and  spec.  The  parameters  are  analogous  as  in  the 
above  charts  for  blochc4  and  spec.  This  program  will  generate  the  spectrum  given  the  parameters 
for  blochc4  and  for  spec.  The  results  are  saved  as  <basefile>x.spec,  <basefile>y.spec,  and 
<basefile>.par.  The  magnetization  components  calculated  for  blochc4  are  deleted  after  the 
spectrum  is  calculated.  This  reduces  disk  space  usage  and  allows  a  more  efficient  way  of 
obtaining  spectra. 

D.  Usage  of  genpar  programs 

The  genpar  type  programs  work  to  generate  a  scripting  file  that  will  perform  multiple  runs  of  the 
bspec  program.  Basic  command  line  is  as  follows. 

Genpar  <Tol>  <linewidth>  <H1>  <T1>  <T2>  <Hm>  <specwidth>  <points> 

<nfreq>  <freqoff>  <basefile> 
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<Tol> 

(real)  Specified  tolerance  of  gaussian  line,  i.e.  number  of  wavepackets  to 
calculate  such  that  on  a  central  line  of  a  gaussian  function,  the  value  is  equal  or 
less  this  tolerance  level. 

<Linewidth> 

Linewidth  of  gaussian  function  (G) 

<H1> 

(real)  Microwave  power  amplitude  (G). 

<T1> 

(real)  Relaxation  Rate  (sec) 

<T2> 

(real)  Relaxation  Rate  (sec) 

<Hm> 

(real)  'A  the  desired  modulation  height  (G)  for  the  simulation.  Okay  to  set  to 
zero  if  no  modulation  is  desired. 

(real)  Width  of  spectrum  (G) 

<Points> 

(integer)  Number  of  points  in  the  spectrum 

('integer')  Number  of  modulation  frequencies  to  calculate 

^ | 

( real)  Modulation  frequency  offset  between  runs 

(filename  string)  Filename  of  the  set  of  generated  spectra. 

E.  Further  Programming  Ideas 

1)  In  the  spec  program,  the  only  linshape  function  (L(H)  as  above)  applied  is  a  gaussian.  More 
functions  would  be  desirable.  Another  improvement  would  be  to  read  in  a  calculated  powder 
spectrum  or  other  generated  spectrums,  and  run  the  spec  program,  interpolating  between  the 
calculated  values  for  the  weight  values  of  each  spin  packet. 

2)  Adaptive  step  sizes  to  reduce  integration  errors  should  be  intoduced.  This  would  allow  a 
more  accurate  simulation  of  the  spinpacket  magnetization  values  as  a  function  of  time,  giving 
a  more  reliable  spectrum  overall.  With  the  current  blochc  program,  the  wavepackets  out  on 
the  wings  are  the  most  inaccurate.  The  inner  wavepackets,  require  a  smaller  step  size  to  be  as 
accurate  as  the  outside,  but  the  constant  stepsize  causes  them  to  overstep.  Using  an  adaptive 
stepsize  would  not  only  increase  the  robustness  of  the  calculations  for  the  outer  packets,  but 
would  most  likely  improve  overall  performance,  due  to  a  smaller  required  step  size  for  the 
inner  packets. 


F.  Source  Listing . 

Bspec2 

#!  /bin/csh  -f 

#1  -  sp  inpackets 
#2  -  subpoints 
#3  -  packet  ponts 
#4  -  Hm 

#5  -  Modulation  Frequency  (Hz) 
#6  -  Microwave  Power  (G) 

#7  -T1  (sec) 

#8  -  T2  (sec) 

#9  -  cycles 
#10  -  spec  points 
#1 1  -  linewidth 
#12  -  begin  field 
#13  -  end  field 
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#14  -  Center 
#15 -Base  File 

,/blochc4  $1  $2  $3  $4  $5  $6  $7  $8  $9  $15 

./spec  $15  $10  $11  $12  $13  $14 
rm  ./work/*  .red 


/*  port  of  Sean  Mcllwain's  matlab  code  for  integrating  the  Bloch  equations. 
* 

* 

*  port  of  file  blochc-m 

*  Sean's  comments  are  preceeded  by  ,m/.f  characters. 

*/ 

/* 

/(Bloch  Equations  for  EPR 

)(liJill  Aproxiamate  the  Bloch  Equations  given  the  parameters 
^Using  Bloch4  Function 
5: 

y.Hl  -  microwave  field 
5: 

'/.  -  Center  of  Resonance 

y.Ti 
y.  te 
*/ 

♦♦include  <string.h> 

♦♦include  <stdio-h> 

♦♦include  <stdlib.h> 

♦♦include  <strings-h> 

♦♦include  <math.h> 

♦♦include  <malloc.h> 

♦♦include  nsutil . hn 
/♦♦♦include  nbdef.hn*/ 

♦♦include  "blochb-h" 


♦♦define  flinPacketStep  O-DS 
♦♦define  FACTOR  IDO 
/♦#def ine  flAXBUF  flO*/ 


/♦  I  suppose  this  is  the  main  program  ♦  / 
void  mainCint  argci  char  ♦argvC]) 

/♦  argvE13  =  spinpackets 
argvEEl  =  substeps 
argvE33  =  filepoints 
argvEMl  =  Hm  Cl/E  Modamp  (Gauss) 
argvCSl  =  v  (Mod  Freq  in  SecA-l) 
argvEt.3  =  HI  (Gauss) 
argvE73  =  T1  (Sec) 
argvCfll  =  TE  (Sec) 
argvE6^  =  cycles 
argv*"103  =  filename 
♦  / 

-C  /♦  beginning  of  main  program  ♦  / 

/* 

/(Parameters  -  Function 
*/ 

long  i*>j; 

long  filepoints; 
long  spinpackets; 
long  points; 
long  subpoints; 
int  Blowup; 

double  modfreq;  /*  /(Hz  */ 

/♦double  pi  =  3 . 14 15TE7 ; */ 

double  Urn;  /♦  ^Modulation  Frequency  (v(Hz)*E*pi)  ♦  / 


double 

Hm ; 

/* 

m/.  Modulation  Amplitude  (Gauss) 

*/ 

double 

ti; 

/* 

/(in 

sec 

*/ 

double 

te; 

/* 

X  in 

sec 

*/ 

double 

hi; 

/* 

/(in 

Gauss  */ 

double 

Ho  =  o; 

/♦ 

‘/.Center  of 

Resonance  (Gauss) 

*/ 

double 

Me  =  i; 

/♦ 

/(Lineshape 

Height  */ 

double 

k  =  i/e; 

/♦ 

‘/.  Lineshape 

Width  */ 

double 

offset ; 

/* 

^Offset 

to  calculate  for  spin 

packet 

*/ 

double 

deltaH ; 

/♦ 

^Number  of  spin 

packets  to 

calculate 

double 

stepsize 

; 

long 

cycles ; 

IV-8 


long  Hc^sizei 
double  *Hci 
double  *timearrayT 
double  maxtimei 
double  ho  =  li 
double  timei 
double  dti 
double  mintimen 
double  temptime1; 
bloch_par  b_p ari 


/* 

‘/.Parameters  -  Integration 
*/ 

/* 

‘/Parameters  -  File  saving 

*/ 

char  param£30D3i 

double  *Hnext i  /*  defined  by  calloc-.  used  as  MCpointslCM]  in  blochM-c  */ 


double  *hsavei 

double  **sresultxi 
double  **sresultyT 
double  **sresultz*i 

double  *tempx i *tempEx t 
double  *tempy *tempBy  n 
double  *tempz n *temp3z i 

/*  assigning  filenames  */ 

char  f iledirtriAXBUFl  =  "./work/":. 

char  f i lenameChAXBUFl t 

char  f i lenameBxth AXBUF3 i 

char  f  i 1 enameByEMAXBUFl i 

char  filename2zEhAXBUF3i 

char  filename3  EhAXBUFli 

FILE  *filelt  *file5i  *file3i 


if  (argc  !=  11) 

■C 

printf  ("Insufficient  #  of  argument s\n" ) i 

printf  ("Usage  blochc  <packets>  <substeps>  <fpts>  <Hm>  <v>"):. 
printf  ("  <H1>  <T1>  <T2>  <cycles>  <f  i lename>\n" ) 
exit (-1) i 

> 

else 

•C 

spinpackets  =  atoi ( argvtll ) n 
subpoints  =  atoi ( argvC23 )  n 
filepoints  =  atoi  (argvC33)i 
points  =  filepointsi 
/* 

Hm  =  atof (argvC4])i 
modfreq  =  at o f ( argvCSH )  i 
HI  =  atof ( argvCbl ) n 
T1  =  atof ( argvC73 ) t 
T2  =  atof ( argvEfil )  i 
*/ 

modfreq  =  a t o f ( argvC53  )  i 

b_par*ho  =  li 

b_p  ar.Hl  =  atof(argvCbl)n 

b_par-Tl  =  atof ( argvE73 ) ^ 

b__par.T2  =  atof  ( argvEfil )  i 

b_par-Hm  =  at of ( arg vC4 3 ) i 

b_par.Um  =  atof ( argvC53 )  *  2-0  *  pii 
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b_par-Ho  =  D • 0 i 

cycles  *  atoi (argvCTl) \ 
strcpyCfilenameiargvClOl)  i 


f ilename2xC03  =  'XD't 

strcpy(filsname2xifiledir ) t 
strcat(filename2xTfilename)T 
strcat(filename2x-i"x-red") t 

f ilename2yE03  =  'NO'i 

strcpyCfilenameSyifiledir)^ 

strcat ( f ilename2y i f ilename) t 

strcat(filename2yTny«redn)  i 

f ilename2zC03  = 

strcpy(filename2zTfiledir ) t 
strcat  ( filen  am  eBz-ifilen  am  e)=i 
strcat(filename2z-i"z*red")T 

filename3C03  =  r\D,:i 

strcpyCfilenameBTfiledir ) i 
strcatCfilename3Tfilename):i 
strcat (f ilename3i" • par") i 

/*  Urn  =  2*pi*modf reqi */ 

maxtime  =  (double) (cycles+1)  *  (1 . 0/modf req ) i 

mintime  =  (l/modfreq)n 

del taH  =  (double)  spinpackets i 

if  ( (l/b_par . T2/HzperGauss )  <  b_par-Hl) 

•C 

stepsize  =  l/b_par - T2/HzperGauss \ 

> 

else 

stepsize  =  b_par.Hli 

> 

if  (MinPacketStep  >=  stepsize) 

i 

stepsize  -  MinPacketStep i 

> 

offset  =  (del taH-1 • 0 ) /2 • 0  *  stepsizei 
Hc_size  =  spinpackets+li 

printf  ( "\nSpinpackets  =  '/.d\nn  t  sp inpackets  ) -. 

printf  ("points  =  Xd\n"  points  )  i 

printf  ("filepoints  =  *d\n"  f  ilepoints  )  ^ 

printf  ("Hm  =  *f  \n"  -»b_par  .  Hm  )  . 

printf  ("liJm  =  */.f  \n"  b_par  .  him )  i 

printf  ("HI  =  Xf  \nn  b_par  .  HI ) 

printf  ("Tl  -  m4f \n" t b__par • T1 ) i 

printf  ( "T2  =  Xf  \n"  i  b_par .  T2  ) 

printf  ("cycles  =  */.d\n" cycles )  i 

printf  ("filename  =  /is\n"  t  f  i  lename  )  **. 

printf  ("maxtime  =  Xf  \n" nmaxt ime ) n 

printf  ("offset  =  #/.f  \n"iof  f set ) 

printf  ("DeltaH  =  Xf  \n" -.deltaH)  i 

printf  ("stepsize  =  Xf  \n" steps  ize  )  i 

printf  ( "Hc_s ize  =  Xd  \n" -.Hc^size )  i 


> 


printf ("ready  to  calloc  \n")i 

/*  from  seanutil*c  */ 
hsave  r  array ( Hc_s ize ) \ 

He  =  array (spinpackets 
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timearray  =  array ( f ilepoints )  i 

sresultx  =  matrix  ( f  ilepoints-.Hc_size>  i 
sresulty  =  matrix ( fi lepoints iHc_size ) i 
sresultz  =  matrix(filepoints-.Hc_size)  i 
dnext  =  array  (4)i 


pr intf ( "done  with  calloc  \n")i 

/*  filling  the  He  array  */ 

for  (i  =  Di  i  <  spinpacketsi  i  =  i  +  1) 

He til  =  b_par-Ho-of f set  +  ((double)  i)  *  stepsizen 
/*printf("Hctkd]  =  kf\nn-iiiHcCi3)i*/ 


> 


/*  assigning  param  string  */ 
/* 

kBuild  parameter  file 


sprint  f(pa'"ai»i"pack.ts\tpoints\tcycles\tHl\tTlNtT3\tmodfreq\tHm\tHo\nXd\tXd\tI'.d\tHgNt29\tH9\tXg\ 

t'i9''cyclesib^partHlib^par  ■  Tl°b!lpar -T2imodf  req'b_p.r-Hmib_par  .Ho)  t 


/*sprintf (param n" 

modfreq=  */.f  Hz\n 
Tl=  Zf  sec\n 
T2  =  ‘/.f  sec  \n 

Hl=  */.f  Gauss  Center=  Zf  Gauss  \n 
Mod  Amp  =  m/.f  Gauss\nn 
imodfreqiTlnTEiHliHonHm/2)  ; 

♦  / 


/♦ 

’/'Calculation  of 

♦  / 

i  =  o; 

Blowup  =  TRUE  t 
time  =  0  ■  0 
flnextCDl  -  0  •  □ ; 
MnextCll  =  □  •  0 ; 
MnextEEl  =  D-0i 
MnextC33  =  Hoi 
j  =  Dn 
print f  ("Starting 


a  spin  packet 


integrat ion\n" ) ; 


while  (Blowup  ==  TRUE  88  i  <  filepoints) 

{ 

Blowup  ”  false; 

dt  =  maxtime  /  (  (filepoints-!)*^  (subpoints) 


while  ((Blowup  ==  FALSE)  &&  (time  <  mintime)) 
I 


/♦ 


♦substep  is  a  procedure  that  affects  time  and  Mnext 

♦  it  basically  calculates  the  next  value(s)  of  the 

♦  magnetization  vector  and  time-,  given  dt*  It  also 
♦has  a  feature  that  it  returns  an  error  value  to 
♦which  signifies  a  blowup  state-  Using  this  states 
♦rest  of  the  program  can  determine. if  it  needs  to 
♦increase  the  #  of  substeps  for  which  to  take 

♦  / 


Blowup  =  substep(f1next-iHcCj]i  b_par  n  &  t  ime  idt-isubpoints); 
>  /♦end  of  while  ((Blowup  ==  FALSE)  &a  (time  <  mintime))*/ 
print f  ("Current  time  is  Z.lg\n"  n  time ) ; 


sresultxEOUCll  =  MnextClB; 
sresultyEQlEll  -  MnextEEJ-. 
/♦sresul tzEDlEl J  =  flnextEBl ;  */ 
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timearrayE03  =  time’-, 


sresul txE03E03  =  timen 
sresultyC03CQ3  =  time1, 

/♦sresultzCG3IC03  =  timei*/ 
i  =  1  i 

printf  ("Entering  Integration  Loop  \n")i 
while  ((Blowup  ==  FALSE)  &&  (i  <  filepoints)) 

/♦printf  ("Filepoint  #  *d  of  Xd\n"  *.  i  *.  filepoints ) \*/ 
/♦printf  ("Calling  substep\n" ) n ♦/ 

Blowup  =  s ubs t ep  ( Unex t -> HcC  j3 1 b_p ar -i a t ime t dt i s ubpo i nts ) -i 

if  (Blowup  ==  FALSE) 

I 

/♦printf  ("time  =  */.f \n"  -i  t  ime )  \  ♦/ 
timearraytil  =  timen 
sresultxfci3E03  =  timen 
sresu 1 t y C i3E03  =  time1, 

/♦sresultzEi3E03  =  timei*/ 

/♦printf  ("sresultxE/.d3E03  =  /Cf\n"iiitime)!i*/ 

/♦if  (ifclOO  ==  0) 

printf  ("Assigning  #’/.d  to  sresultE’/dl\n" •,  i  n  j  +  1 )  \ 

y  ♦/ 

sresultxCi3C13  =  HnextElli 
sresultyE  ilEll  =  (InextCEln 
/♦sresultzE i3C13  =  f1nextE33i*/ 
i  +  +  t 


>/♦  end  of  while  ((Blowup  ==  FALSE)  aa  (i  <  f i lepoints ) ) ♦/ 

if  (Blowup  ==  TRUE) 

-C 

i  =  o* 

time  =  0  •  0 *-> 
f1nextE03  =  O.On 
f1nextE13  =  0-0i 
MnextC23  =  O-On 
f1nextE33  =  hoi 

j  =  o; 

printf  ( "\nBlew  up-,  trying  V.6  +  *d"nSubpoints  -.FACTOR )  n 
subpoints  =  subpoints  +  FACTORS 
printf  ("  -  /.d\n"-,subpoints):. 

> 

> 

//Now  have  a  substep  that  doesn't  cause  the  program  to  blow  up 
//generate  cos  array  for  mass  scale  integration. 

//gencos  (b_par  .  UJm  ->dt )  \ 

printf  ("Starting  second  integration^")  n 

for  ( j=li j<sp inpackets i j  +  +  ) 

■C 


time 

=  O-On 

flnex  tE03 

=  On 

MnextE13 

=  On 

flnex  tE23 

=  o; 

MnextE33 

=  b_par  •  flo  n 

/♦ 

♦calculate  M  up  to  min  time  in  hopes  of  removing  the 
♦transient  from  the  spectrum 
♦  / 

printf  ("Advancing  time  to  ‘/g\n"imintime)i 
while  (time  <  mintime) 

/♦ 


IV-12 


*substep  is  a  procedure  that  affects  time  and  Mnext 
♦  it  basically  calculates  the  next  value (s )  of  the 
*magnet ization  vector  and  timei  given  dt-  It  also 
*has  a  feature  that  it  returns  an  error  value  to 
*which  signifies  a  blowup  state*  Using  this  statei 
*rest  of  the  program  can  determine  if  it  needs  to 
♦increase  the  #  of  substeps  for  which  to  take 
*/ 

Blowup  =  subst ep  ( Unex 1 1  HcC  j3  t  b_par  t  & t  ime  i  dt  i subpo i nts  )  •» 

>  /*end  of  while  ((time  <  mintime))*/ 
printf  ("Current  time  is  *lg\n"  *. time )  \ 
sresultxED3E j+13  =  HnextE13n 
sresultyE03E j+13  =  MnextEE3i 
sresfjltzE03Ej  +  13  =  MnextE33n 


/♦printf  ("Entering  Integration  Loop  Nn")1.*/ 


printf  ("Into  integration^" )  i 
cosarray  =  array  (subpoints+1)  i 
cosarrayh  =  array  ( subpo int s  )  t 

for  (i=lii  <  pointsni++) 

I 

//printf  ("point  #  kd\n"ii)i 

//printf  ("time  ’/.lg\n" i t imearrayE i'3 )  i 

temptime  =  t imearrayE i3 i 

tempx  =  sresul txEi-13  i 

tempy  =  sresultyEi-13i 

tempz  =  sresultzEi-13:i 

tempEx  =  sresul  txEil  =. 
tempEy  =  sresul tyE il  i 
tempEz  =  sresultzEili 

gencosE  ( mod f  req  b_par  .  Hm  i  tempt  i me  dt  subpoint s  +  1  cosarray ) 
gencosE  (modfreqi b_par .Hm-itemptime  +  0-  5*dtndtTSubpoints-iCosarrayh) 

for  (j  =  li  j  <  spinpacketsi j++) 

time  =  temptime^ 
f1rextED3  =  timei 

//MnextE13  =  sresul txEi-13E j+13 i 
//MnextEEI  =  sresul tyEi-13E j  +  13i 
//flnext E33  =  sresul tzEi-13E j+13 ^ 

f1nextE13  =  tempxEj+13i 
HnextEE3  =  tempyEj+13i 
MnextE33  =  tempzEj+13i 

/♦printf  ("Filepoint  #  of  5cd\n"  n  i  t  f  i lepoints )  t  ♦/ 

//printf  ("Calling  substep\n")^ 

Blowup  =  substep(nnext-.HcEj3-.b_par-.&time-.dt’iSubpoints)*-, 

/♦if  (i/.100  ==  D) 

printf  ("Assigning  #/.d  to  sresultE Xd3\n" ■»  i  «.  j  +  1 )  ^ 

>  ♦/ 

//printf  ("Assigning  values\n")i 
//sresultxEi3E  j  +  13  =  MnextE13:i 
//sresultyEi3C j+13  =  MnextEE3i 
//sresul  tzEi3E  j+13  =  MnextE33i 
tempExE  j  +  13  =  hnextC13i 
tempEyEj+13  -  HnextEE3n 
tempSzE j+13  =  MnextE33i 

>/♦  end  of  for  (j  <  spinpackets ) ♦/ 

>  /♦end  of  for  (i  <  points)*/ 
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hsaveCDl  =  Di 

for  (i  =  li  i  <  (Hc_size);  i  =  i  +  1) 

hsavetil  =  Ho-HcEi-132.  /*  fcPlace  Delta  as  a  header  */ 

> 


/* 

V.Save  Files  ^  4 

/Ceval  (['save  f  filenamelx  1  savefilelx  -ascii  -tabs  -double  j)t 
*/ 

printf(ndeltaH  =  %f  \nnideltaH)T 
printf  ("f  ilepoints  =  */.d  \nn  .  f  ilepoints )  i 


printf(nnow  opening  filel  \n")i 
f i lei  =  fopen(f  ilenameEx-."w")  n 
if  (Ifilel) 

printf  ("filel  opening  seems  to  have  failedXn")1! 
exit (-1) n 

> 

else 

printf ("filel  opening  seems  to  have  succeeded  \n 

for  (j  =  Oi  j  <  Hc_size  =i  j  +  +) 

-C 

fprintf  (  f  i  lei -i  "*g\t" hsavet  jl)  ^ 

> 

fprintf(fileln"\n")*T 

for  (j  s  Di  j  <  filepointST  j  =  j  +  1) 


/*printf("j  =  Zd  \n"-.j)i  /*  */ 


/*fprintf(filel-i"*g\t"-isresultxCj3C(int)deltaH])n*/ 


for  (i  -  Di  i  <  H c_s izei  i  =  i  +  1) 

-C 

/*  printf("In  Loop\nn) i 

printf ("sresult\i*d3E*d3  =  ‘/.f  \n" j  ■.  i  *. sresul txE  j3Ei3 ) 

*/ 

fprintf  (fileli"*g\t"isresultxEj3Ei3)  t 

> 

fprintf  ( filel i"\n") ^ 

> 

f close (filel)  n 


/Ceval  (E'save  1  filenamely  1  savefilely  -ascii  -tabs  -doublef3)i 
*/ 

printf("now  opening  fileE  \n")n 
fileE  =  fopen(f ilenameEyi"w")  i 

/*  fprintf  ( fileS-."*f  *f \n" n hsaveE03 i hsaveE13 )  n  /*  */ 
for  (j  -  Oi  j  <  Hc_size  \  j++) 

•C 

fprintf  (fileE*i")Cg\t"ThsaveEj3)  t 

> 

fprintf(fileli"\n")  i 

for  (j  =  Di  j  <  filepointSn  j  =  j  +  1 ) 

/*fprintf(fileE-i"^g\tnTsresultyEj3E(int)deltaH3)T*/ 

for  (i  =  Oi  i  <  Hc^size;  i  =  i  +  1) 

■C 

f printf (fileE, "Xg\t" sresul tyE j3Ei3)  i 
> 

fprintf(fileSi"\nn) i 
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> 

f  close  ( f ileE )  =. 


/Ceval  (C'save  '  filenamelz  1  savefilelz  -ascii  -tabs  -double'!) 
*/ 

/* 

printf("now  opening  file3  \n")i 
file3  =  fopenCf  i  lenameEz "w" )  i 

for  (j  =  Oi  j  <  Hc_size  =.  j  +  +) 

-C 

fprintf(file3in/Cg\tnThsaveCj!)T 

} 

fprintf(fileli"\nn)i 

for  (j  =  Di  j  <  filepointsn  j  *  j  +  1) 

■C 

for  Ci  =  Qi  i  <  Hc_sizei  i  =  i  +  1) 

fprintf  (f  ile3in>'.g\tn-1sresultzC  j!Ci!)  i 

> 

fprintf(file3in\nn)i 

> 

f close ( f i le3 )  t 
*/ 

printf("now  opening  filel  again  \n")i 
filel  =  fopenCf ilename3nnwn) t 
fprintf(fileli"/'.s\nnTparam)i 


f close ( f i lei )  t 

f ree_matr ix  (sresultxnfilepoints) ^ 
f ree_matrix  (sresultyTfilepoints)n 
/*  f r ee_ma t r i x  (sresultzifilepoints)!*/ 
f r e e_a rray  ( fl  n  e  x  t )  n 

>  /*  end  of  main  program  */ 
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